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CHAPTER  1   INTRODUCTION 

It  is  a  complex  and  uncertain  world.  Probability  is  used  everywhere 
today.  But,  the  word  "probability"  is  often  accepted  rather  skeptically 
by  engineers;  they  think  that  they  are  employing  an  "exact  science"  despite 
much  use  of  empiricism.  We  can  say  that  engineering  is  a  type  of  applied 
physics.  Many  physical  phenomena  can  be  described  only  by  using  the  theory 
of  probability.  It  is  not  surprising  to  discover  that  equipment  failures 
which  result  from  interactions  of  heart,  electric  and  magnetic  field, 
"static"  loadings,  and  vibrations  can  best  be  described  in  probabilistic 
terms  [16]. 

Some  of  the  engineering  design  are  still  stuck  on  conventional  design 
methodology  by  using  "safety  factors".  It  should  be  brought  to  the  area 
of  probabilistic  design.  Reliability  is  a  terminology  in  the  area  of 
probabilistic  design. 

What  is  rel  i ability?  There  are  several  definitions  of  reliability; 
one  of  the  best  is  that  of  the  National  Aeronautics  and  Space  Adminis- 
tration (NASA)  [2]  which  defines  "reliability  "as"  the  probability  of  a 
device  performing  adequately  for  the  period  of  time  intended  under  the 
operating  conditions  encounted."  Other  definitions  are:  "reliability 
is  a  meausre  of  the  capacity  of  a  piece  of  equipment  to  operate  without 
failure  when  put  into  service"  [16],  and  "The  reliability  of  a  system 
is  the  probability  that,  when  operating  under  stated  environmental  con- 
ditions, the  system  will  perform  its  intended  function  adequately  for 
a  specified  interval  of  time"  [10]. 


In  this  thesis,  reliability  is  defined  as:  the  probability  of 
the  strength  of  a  component,  which  is  a  certain  stress  resisting  capacity 
of  the  component,  exceeding  the  stress  induced  by  the  operating  con- 
ditions. If  the  stress  exceeds  this  capacity,  then  failure  results.  It 
is  unfortunate  to  use  words  such  as  "stress"  or  "strength".  The  imputed 
meaning  is  not  necessarily  restricted  to  mechanical  loading  in  the  minds 
of  engineers.  We  can  use  it  in  a  broader  sense,  applicable  in  a  variety 
of  situations  well  beyond  the  traditional  mechanical  or  structural 
systems.  "Stress"  is  used  to  denote  any  agency  that  tends  to  induce 
"failure",  while  "strength"  denotes  any  agency  that  tends  to  resist 
"failure".  "Failure"  itself  is  taken  to  mean  failure  to  function  as 
intended;  it  is  defined  to  have  occured  when  the  actual  stress  exceeds 
the  actual  strength  in  the  beginning  [10]. 

Mathematical  formulation  of  reliability  is  given  by 

Reliability  =  R  =  P(S  >  s)  =  P(S  -  s  >  0)  (1.1) 

where  S  is  strength,  s,  stress,  and  P,  probability  function. 

The  detailed  expression  of  reliability  will  be  discussed  in  Chapter  2. 

Literature  survey:  This  thesis  is  mainly  concerning  about  optimizing 
the  following  objectives;  maximize  reliability  of  stress-strength  inter- 
ference model,  minimize  total  cost  of  a  component,  minimize  weight  of  a 
component;  subject  to  some  absolute  constraints.  The  problem  will  be 
solved  by  multiple  objective  decision  making  methods.  The  idea  of  this 
thesis  is  developed  according  to  the  following  literature  survey. 

Since  1961,  I.  Bazovsky  [2]  has  mentioned  the  stress-strength 
interference  theory  (S-S-I).  In  1964,  M.  J.  Bratt,  G.  Reethof,  and 


G.  W.  Weber  [3]  interpreted  the  stress-strength  interference  case  and 
also  developed  some  models  for  time  varying  stress-strength  interference 
cases  (S-S-T).   In  1964,  R.  L.  Disney,  C.  Lipson,  and  N.  J.  Sheth  [4] 
developed  the  complete  models  of  S-S-I  of  several  different  distributions 
In  1972,  Leonard  Shaw,  Martin  Shooman,  Robert  Schutz  [15]  developed  the 
several  models  for  S-S-T  models.   In  1975,  K.  C.  Kupar  [9]  studied 
the  optimization  technique  of  a  component  reliability  problem  witn  a 
cost  constraint.  The  problem  dealt  mainly  with  one  objective,  sub- 
ject to  one  or  more  constraints.  But  in  the  real  world,  the  designer 
or  decision  maker  would  like  to  optimize  several  objectives  at  the  same 
time  rather  than  one  objective. 

Objective:  The  objective  of  this  thesis  is  to  solve  the  following 
multiple  objectives  problem.  To  find  x_  so  as  to 
Max  fR(x) 

Min  fc(x) 

Min  fw(x) 

subject  to   g_(x)  <_   0 

where  fR,  fp,  and  f,.  are  the  functions  of  reliability,  cost,  and  weight 
respectively,  x  denotes  a  decision  variable  vector  (x-,,  x?,  ...,  x  ). 

In  solving  this  problem,  it  is  necessary  to  replace  the  concept  of 
optimum  with  that  of  best  compromise.   The  problem  will  be  solved  by 
two  methods  of  the  Multiple  Objective  Decision  Making-methods  (MODM) 
[6],  goal  programming  [7,9]  and  sequential  multiobjective  problem 
solving  technique  (SEMOPS)  [14]. 


In  order  to  develop  the  model  as  we  mentioned  above,  Chaper  2  to 
Chapter  4  provide  the  background  for  getting  the  model.  Chapter  2  dis- 
cusses the  distributions  of  stress  or  strength,  the  relationship  between 
safety  factor  and  reliability,  and  computation  of  moments  of  a  function  of 
random  variables.  Chapter  3  discusses  several  stress-strength  inter- 
ference models.  Chapter  4  discusses  several  time  dependent  stress- 
strength  interference  models.  In  Chapter  5,  the  main  part  of  this  study, 
development  of  the  problem,  and  how  to  use  different  proper  methods  to 
solve  the  problem  are  presented.  Chapter  6  gives  the  conclusion  and  the 
extensions  of  the  problem  which  remain  unsolved. 


CHAPTER  2.  PROBABILISTIC  ENGINEERING  DESIGN 

The  conventional  design  approach,  which  is  based  on  safety  factors 
and  safety  margins,  gives  little  indication  of  the  probability  of  failure 
of  the  component.  Some  designers  believed  that  a  component  failure  could 
be  completely  eliminated  by  using  a  safety  factor  above  a  certain  large 
magnitude.  Actually,  the  failure  probability  may  vary  from  a  low  to  an 
extremely  high  value  for  the  same  safety  factor.  The  safety  factor  can 
be  used  only  when  its  value  is  based  on  considerable  experience  with 
parts  similiar  to  the  one  considered.  Furthermore,  the  design  variables 
and  parameters  are  often  random  variables,  which  are  ignored  by  the  con- 
ventional design  approach.  Therefore,  another  design  methodology  which 
does  consider  the  probability  nature  of  the  design  is  needed.  Such  a 
design  methodology  is  called  "probabilistic  design."  It  identifies 
explicitly  all  the  design  variables  and  parameters  which  can  determine 
both  the  stress  and  strength  distributions  (Fig.  1).  When  these  two 
distributions  are  determined,  the  component  reliability  can  be  calculated 

2.1  Distributions  of  stress  and  strength 

For  the  strength  computations,  considerations  must  be  given  to 
the  properties  of  the  material  used.  The  factors  which  affect  the 
strength  may  include,  but  are  not  necessarily  limited  to,  the  following: 
(Fig.  1)  [12]. 
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Fig.  1  Stress  and  strength  factors  and  distributions  [11] 
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The  factors  which  affect  the  stress  may  include,  but  are  not 
necessarily  limited  to,  the  following:  (Fig.  1)  [12] 

1.  Stress  concentration  factors 

2.  Load  factors 

3.  Temperature  stress  factors 

4.  Forming  or  manufacturing  stress  factors 

5.  Surface  treatment  stress  factors 

6.  Heat  treatment  stress  factors 

7.  Assembly  stress  factors 

8.  Corrosion  stress  factors 

9.  Direct  surface  environment  stress  factors 
10.  Notch  sensitivity  factors. 

Here,  the  determination  of  stress  and  strength  should  be  limited  to 
mechanical  components  where  the  fiber  stress,  in  pounds  per  square  inch, 
is  the  stress  that  governs  failure.  The  theory  that  is  developed  is  not 
limited  in  application  to  this  case  only. 

The  distributions  of  stress  and  strength  generally  can  be  classified 
into  the  following  distributions: 


(1)  Normal  distribution:   (Fig.  2) 

The  probability  density  function  of  a  normal  distribution  can  be 
given  by 

-  2> 


where 


f(x)  =  — —  exp  -  j  (— ) 


x  =  mean  value  of  variable  x 


a  =  standard  deviation  of  variable  x 
x 


00   <   X   <   °° 


(2.1) 


For  example,  the  probability  distribution  of  the  length  of  human  lives 
is  found  to  be  normal.  Also,  the  probability  distribution  of  ultimate  ten- 
sile, yield  and  endurance  strengths  of  steels  are  found  to  be  normal.  The 
probability  distribution  of  the  stress  such  as  rocket  motor  thrust  or  the 
gas  pressure  in  the  cylinder  heads  of  reciprocating  engines  are  also  found 
to  be  normal . 

(2)  Log-normal  distribution:  (Fig.  3) 
The  probability  density  function  is  given  by 


f(x)  =  —  exp 2  (£n  x  '  v)    \    ,       x  >  0 

Xa*'2T7     I   2o  J 


(2.2) 


where  y  and  o  are  the  mean  and  the  standard  deviation  respectively  of  the 
variable  zn   x,  which  is  normally  distributed.  If  we  let  y  =  £n  x,  then 

dy  =  (— )dx.  From  equation  (2.2),  we  have 


f(y)  =  -I- exp  -  \   (y  -  y)2 

a/2-rr     *•   2a  -* 


y 


co  <  V  <  °° 


(2.3) 


It  is  found  that  the  strength  properties  of  structural  alloy  materials 
often  tend  to  follow  a  log  normal  distribution. 


Figure  2     Normal    distribution 
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Figure  3.     Log  normal   distribution  when  a1  =  1 
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(3)  Exponential  distribution:   (Fig.  4) 
The  probability  density  function  is  given  by 

f(x)  =  Xx  e  x  ,       0  <  x  <  co  (2.4) 

where  the  mean  value  of  x,  p,  and  the  standard  deviation  of  x,  a,  was  the 
following  relations: 

1       .   1 

y  =  r  '    G "  r 

X  X 

(4)  Gamma  distribution:   (Fig.  5) 

The  probability  density  function   is  given  by 

,n     n-1     -}c: 

fM   =         *w7TT •         n  >  °>       *  >  °>       0  <  x  <  -  (2.5) 

I  [n )  — 

where  a  is  called  the  scale  parameter  and  n  the  shape  parameter  and  r(n)  is 
the  gamma  function. 

(5)  Weibull  distribution 

The  probability  density  function  for  a  Weibull  distribution  x  is  given  by 

R     x-xo  6-1     r       x-xo  5) 
f(x)  =  (f)(-r^)     exp  -  (— &)    x  ■:  x  <  -         (2.6) 


e 


o  - 


where   6  is  the  slope  parameter,  6  is  the  scale  parameter,  and  x  is 

the  location  parameter  or  lower  bound  of  x. 
The  distribution  becomes  the  exponential  when  a  =  1,  a  distribution 
skewed  to  the  left  when  1  <  a  <  3,  the  normal  when  a  =  3,  approximately, 
and  a  distribution  skewed  to  the  right  when  a  >  3. 
(6)  The  extreme  value  distributions: 

In  general,  the  extreme  value  distributions  are  applicable  where  the 
phenomena  causing  failure  depend  on  the  smallest  or  the  largest  value  from 
a  sequence  of  random  variables. 


12 


Figure  4.  Exponential  distribution 
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Figure  b     Gamma  distribution 
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There  are  two  type  of  distributions  of  extreme  value,  one  is  the 

smallest  value  distribution,  another  is  the  largest  value  distribution. 

x 
If  the  cumulative  function  associated  with  f(x)  is  F(x)  =  /  f(y)dy,  then 

—00 

the  distribution  of  the  smallest  value,  y  =  min(x-,,x?,  ...,  x  ),  in  the 
sample  size  n  drawn  from  the  population  is  given  by  the  probability  density 
function 

gniy)  =  nf(y)(l-FCy))n'1  (2.7) 

or  by  the  cumulative  distribution  function 

oo 

Gn(y)  =  /  gn(y)dy  =  1  -  (1  -  F(y))n  (2.8) 

—  00 

Also,  the  distribution  of  the  largest  value,  z  =  max  (x-, ,  x?,  ...,  xp), 
in  the  sample  size  n  drawn  from  the  population  is  given  by  the  probability 
density  function 

gn(z)  =  nf(z)(F(z))n"]  (2#g) 

or  by  the  cumulative  distribution  function 

Gn(z)  =  f   gn(z)dz  =  (F(z))n  (2.10) 

—  oo 

The  distribution  of  fatigue  of  metals  or  textiles,  the  breakdown  of 
dielectrics  and  the  corrosions  of  metals  tend  to  follow  extreme  value 
distribution. 

2.2  General  expression  for  component  reliability  [10]. 

Let  f  (•)  and  fA-)   be  the  density  functions  of  the  stress  (s)  and 
the  strength  (S)  respectively,  (see  Fig.  6.)  Then,  by  definition, 
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Figure  6     Stress-strength  interference 
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Reliability  =  R  =  P(S  >  s)  =  P(S  -  s  >  0)  (2.11) 

The  shaded  portion  in  Fig.  6  is  the  interference  area,  which  is  the  in- 
dication of  the  probability  of  failure.  Let  us  magnify  this  interference 
area  as  shown  in  Fig.  7. 

The  probability  of  a  stress  value  lying  in  a  small  interval  ds  is 
equal  to  the  area  of  the  element  ds  under  the  curve,  that  is 

P(sQ  -  ^f  <  s  <  sQ  +  4|)  =  fs(sQ)ds  (2.12) 

The  probability  that  the  strength  S  is  greater  than  a  certain  stress  s 
is 

oo 

PCS  >  sQ)  =  /  fs(S)dS  (2.13) 

o 
It  is  assumed  that  the  stress  and  the  strength  random  variables 

are  independent  and  greater  than  or  equal  to  zero.  The  probability  for 

the  stress  value  lying  in  the  small  interval  ds  and  the  strength  S  exceeding 

the  stress  given  by  this  small  interval  ds  is  given  by 

00 

fs(s  )  /  fs(S)dS  (2.14) 

s 

0 

Because  event  A  and  event  B  are  independent,  then 

P(AAB)  =  P(A)  •  P(B)  (2.15) 

Now,   the  reliability  of  the  component  is  the  probability  that  the 
strength  S  is  greater  than  the  stress  s  for  all   possible  values  of  the 
stress  s  and  is  given  by 


—  00 


s-'l   i  -S 


R  =  /     Us)     /  fc(S)dS|ds  (2.16) 
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Figure  7.  Enlarged  portion  of  this  interference  area 
of  Figure  6. 


Similiarly,  the  reliability  of  a  component  can  be  expressed  by 
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R  =  /  -^s Cs ) !  J  f  (s)ds 

CO  V   —00 


dS 


(2.17) 


The  unreliability,  R,  is  defined  as 


R  =  probability  of  failure  =  l-R=P(S<_s) 
Substituting  for  R  from  equation  (2.16)  yield 


(2.18) 


co         r       oo  \ 

R  =  P(s  <  s)  =  1  -  /  f  (s)  /  fs(S)dS  ds 

—  oo  ^.   C  J 


=  1  -  /  fs(s)[l  -  F$(s)]ds 


=  /  Fs(s)  fs(s)ds 


(2.19) 


Similiarly,  using  equation  (2.17)  we  have 


R  =  P(S  <  s)  =  1  -  /  fs(S)|/  fs(s)ds]dS 

—  CC  '  ^00 

oo 

=  1  -  /  fQ(S)  •  F(S)dS 


=  /  [1  -  F  (S)]  fs(S)dS 


(2.20) 


Let  y  =  S-s,  then  y  is  called  the  interference  random  variable.  Now 
we  can  define  the  reliability  as  R  =  P(y  >  0).  The  density  function, 
f  (• ) ,  is  then  given  by 


19 


fv(y)  =  /  fc(y  +  s)f  (s)ds  (2.2i) 

y     s        s 

Hence  the  probability  of  failure,  which  is  y  <_  0,  is  given  by 

0         0  « 
R  =  /  f  (y)dy  =  /  /  fs^+s)  .  fs(s)dsdy  (2.22) 

-co   "^  — oo  -y 

and  the  reliability  by,  which  is  y  >  0,  is  given  by 


R  =  /  f  (y)dy  =  /  /  f.Cy+s)  •  f .(s)dsdy  (2.23) 

0  y       0  0   5 

2.3  Definition  of  safety  factor  and  safety  margin  [11] 

A  variety  of  definitions  for  safety  factor  has  been  proposed  and 
used.  These  definitions  are  given  as  follows: 

1.  Ratio  of  ultimate  strength  in  a  component  to  the  allowable  or  actual 

working  stress.  This  can  be  written  as 

S 
Safety  factor  =  F  =  -&  (2.24) 

sw 

where    S  =  ultimate  strength, psi 

s,,  =  workina  stress,   psi 
w 

2.  Ratio  of  yeild  strength  in  a  component  to  the  allowable  or 
actual  working  stress. 

3.  Ratio  of  maximum  safe  load  to  normal  service  load. 

4.  Safety  factor  is  the  ratio  of  computed  strength  S,  to  the  corres- 
ponding computed  load  l. 

5.  Ratio  of  damaging  stress  to  maximum  known  working  stress. 

6.  Ratio  of  mean  strength  to  mean  load. 

7.  Ratio  of  significant  strength  to  significant  stress. 
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It  may  be  seen  that  although  definitions  vary  widely  with  appro- 
priate interpretations  of  the  terms  used,  the  safety  factor  is  the  ratio 
of  a  particular  strength  value  to  a  particular  stress  value.  The  central 
safety  factor  is  the  ratio  of  the  central  tendency  measure  of  the 
strength  distribution  to  the  central  tendency  measure  of  the  stress  dis- 
tribution. 

The  best  concept  such  a  safety  factor  can  convey  is  how  far  the 
mean  strength  is  removed  from  the  mean  stress.  Sometimes  designers 
believe  that  designing  to  a  safety  factor  above  some  preconceived  mag-' 
nitude,  usually  above  2.5  would  result  in  no  component  failure.  On  the 
contrary,  with  such,  and  even  higher  safety  factors,  the  failure  prob- 
ability may  vary  from  a  very  satisfactory  low  value  to  an  intolerably 
high  value.  Furthermore,  a  safety  factor  of  one,  to  most  designers, 
implies  that  failure  will  occur  100%  of  the  time  because,  presumably, 
there  is  no  safety  margin,  whereas  failure  would  actually  occur  only 
50/c  of  the  time,  if  the  stress  and  strength  distributions  are  normal. 

The  basic  safety  margin  definition  used  in  most  design  books  is 

Safety  Margin  =  M  =  F  -  1  (2.25) 

This  means  that  the  safety  margin  is  the  amount  by  which  the  safety 
factor  value  exceeds  unit.  This  definition  again  ignores  the  fact 
that  stress  and  strength  are  distributed,  as  does  the  safety  factor. 

2.4  Safety  factor  and  reliability 

From  Table  1,  we  can  see  that  the  same  safety  factor,  2.5,  can  re- 
sult the  reliability  from  0.987  to  1.0  and  even  the  safety  factor,  5.0, 
may  only  result  in  a  reliability  of  0.9738.  It  is  not  necessarily  true 
that  the  higher  value  of  safety  factor  can  result  in  a  higher  reliability. 
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Table  1     Safety  factors  and  reliability*     [11] 


CASE 

No. 

MEAN 
STRENGTH 

^S 

MEAN 

STRESS 

Ms 

STRENGTH 
STANDARD 
DEVIATION 

Gs 

STRESS 

STANDARD 

DEVIATION 

°s 

FACTOR 

OF 

SAFETY 

RELIABILITY 
R 

1 

50,000 

20,000 

2,000 

2,500 

2.5 

1.0 

2 

50,000 

20,000 

8,000 

3,000 

2.5 

0.9997 

3 

50,000 

20,000 

10,000 

3,000 

2.5 

0.9979 

4 

50,000 

20,000 

8,000 

7,500 

2.5 

0.9965   " 

5 

50,000 

20,000 

12,000 

6,000 

2.5 

0.987 

6 

25,000 

10,000 

2,000 

2,500 

2.5 

0.9(6)4 

7 

25,000 

10,000 

1,000 

1,500 

2.5 

oJ]6h 

8 

50,000 

10,000 

20,000 

5,000 

5.0 

0.9738 

9 

50,000 

40,000 

2,000 

2,500 

1.25 

0.99909 

10 

50,000 

10,000 

5,000 

5,000 

5.0 

1.0 

Assume  that  the  stress  and  strength  are  normally  distributed. 


22 


The  relationship  between  safety  factor  F  and  reliability  will  be  given 
as  follows: 

Letting  safety  factor  F  be  a  random  variable,  we  define  F  as  the  ratio 

S       - 
of  the  strength  S  to  the  stress  s; that  is, F  =  -  .  Let  S  and  s  denote 

the  mean  strength  and  stress,  respectively.  Let  F  denote  the  central 
factorof  safety,  which  is,  F  =  S/s.  Using  the  Bienayme-Chebyshev  in- 
equality and  Taylor's  series  approximation,  we  can  get  the  following  re- 
lationship between  central  safety  factor  F  and  reliability  R  [10]. 


Fci 


where 

cs  °S 

vs  =  r>Vr  (2-26) 

2.5  Computation  of  moments  of  a  function  of  random  variables 

Let  us  consider  the  case  where  x  is  a  one-dimensional  variable. 
The  expression  of  Y  =  f(x)  about  the  point  x  =  y,  which  y  is  expected 
value  of  variable,  x,  by  Taylor's  series  approximation  up  to  the  first 
three  terms  is 

t        \2 
Y  =  f(x)  =  f(y)  +  (x-y)f'(y)  +  &ff-   f"(y)  +  r  (2.27) 

where  r  is  the  remainder.  Taking  the  expectation  of  equation  (2.27), 
we  have 
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E(Y)  =  E[f(v)]  +  E{xf'(v)  -  yf'(y)) 
+  E{1  fn(v)(x-v#  +  E(r) 
=  f(y)  +  (yf'(y)  -  uf(v))  +  \  fH(v)V(x)  +  E(r) 

=  f(v)  +  |f"(y)V(x)  (2.28) 

Equation  (2.28)  is  an  approximation  for  the  expected  value  of  Y 
because  we  have  ignored  the  remainded  terms  in  the  Taylor's  series  ex-, 
pansion.  If  the  variance  of  x  is  small,  then  we  may  further  ignore  the 
second  term  in  equation  (2.28)  to  obtain 

E(Y)  =  E[f(x)]  ~-   f(y)  (2.29) 

The  approximation  value  of  V(Y),  is  obtained  by  considering  once 
again  the  Taylor's  series  expansion  up  to  the  first  two  terms: 

Y  =  f(y)  +  (x-y)f'(y)  +  r  (2.30) 

Taking  the  variance  of  equation  (2.30)  and  ignoring  the  remainder  term  r, 
we  have 

V(Y)  =  V[f(v)]  +  VECx-n)f'(p)] 

-  [f (y)]2V(x)  (2.31) 

If  we  consider  x  is  n-dimensional   variables,   that  is 

Y  =  f(x-,,  x2.    ...,  xn)   =  f(x)  (2.32) 

Let  y  =  (y-,,  •-.,  y  )  and  o  =  (a,,  ...,  c   )  denote  the  vectors  of 
the  expected  values  and  the  standard  deviations  of  x-, ,  x^,  •••,  x 
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respectively.  Then,  by  the  Taylor's  series  expansion,  we  have 


Y  =  f(x-j,  x2,  ...  xn) 


=  f(y-,  ,y?,  ...  y  )  +  I 


n  3f(x) 


1^2 


(x,-u,) 


n    >-|  3xi   'x=y 


.   n   n  3  f(x)  I 

+  2T  J,  ,L  jr^r  lx-u  (v^'W  +  r 


j=l  i=l  0AiCAj 


(2.33) 


Taking  the  expectation  of  equation   (2.33),  we  have 


n     3f(x)    I 
E(Y)  =  f(Vl,    ....pn)+     l    -^-  \x__v  EU^) 


2 
,      n       n     3   f(x) 


x=yE[(xry.)(xry.)]  +  E(r) 


(2.34) 


If  X-, ,  x«s  ...,  x  are  independent  variables,  then,  after  deleting 
the  zero  terms,  we  get 


E(Y)  =  f(yr  ...  yn)  +  £  _£ 


1  n  s  f(x) 


3x?    X=y    ' 


V(x.)  +  E(r) 


,  n  3^f(x) 
-  f(vi.  ...  yJ  +  i  A 


l 


V    2i=l   3X^ 


V(x.) 

X  =  p   x  1  ' 


(2.35) 
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If  we  ignore  the  second  term  in  equation  (2.35)  then, 


E(Y)  =  f(yr  ...  un) 


(2.36) 


Now,  considering  only  the  first  two  terms  in  equation  (2.33)  and  taking 
the  variance,  we  have 


n  Sf(x) 

V(Y)  =  V[f(x)|   ]  +  V|  I 

x-y     li  =  1  ox. 


r-v   txrp,) 


n 

=  y 

i=l 


r 


of(x) 


12 


x=y 


V(x.) 


(2.37) 
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CHAPTER  3     THE  STRESS  AND  STRENGTH   INTERFERENCE  MODELS   (S-S-I) 

In  order  to  compute  the  reliability  we  have  to  know  the  nature  of  the 
stress (s)   and  strength(S)   random  variables.     As  we  mentioned  before  in 
Section  2-1,  the  distributions  of  stress  or  strength  can  be  a  normal, 
log-normal,   exponential,  gamma,  Weibull,  or  extreme  value   distribution. 

We  can  use  equations   (2.16)  and   (2.17)   to  compute  the  reliability. 
Tney  are 

00  /•       00  \ 

R  =  /     fs(s)     /     fs(S)dSjds 


s 


or 


R  =  /     fs(S)[/     fs(s)as]dS 

—  00  *.  —  00  J 

where  f  (•)  and  f<.(.)  are  the  prooability  density  functions   (pdf)  of 
stress  and  strength,  respectively.      They  can  be  one  of  distributions  we 
mentioned  above.     Hence,   the  reliability  of  a  component  can  be  a  dif- 
ferent type  of  formulation.     Tnere  are  many  combinations  of  different 
distributions.     But,  only  some  of  them  are  derived,   and  the  rest  of 
them  can  be  seen  in  references   [4]  and   [10]. 

3.1     Reliability  computation  for  normally  distributed  strength  and 
stress  [10]. 
The  pdf  of  a  normally  distributed  stress(s)   is  given  by 

f  (s)  =— L-exp   I-  \\r^\  -  -  <  s  <  -  (3.1) 
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fc(S)  -  — *--  exp  I-  i  [— A]  -  .  <  S  <  -     (3.2) 


and  the  pdf  of  a  normally  distributed  strength  (S)  is  given  by 

2 
S_- 

where 

us  =  mean  value  of  strength 
m  =  mean  value  of  stress 

a~  =   standard  deviation  of  strength 
a     =  standard  deviation  of  stress 

Let  us  define  y  =  S-s.  The  random  variable  y  is  also  normally  dis- 
tributed with  mean  y  ,  that  is 

^y  =  US  '  ps  (3'3) 

and  a  standard  deviation  o  is 

y 


ay  =]    os2  +  cs2  (3.4) 


The  reliability  R  can  be  expressed  in  terms  of  y  which  y  >  0,   by 
R  =  P(y  >  0) 

'  /    —^exp   !-  l[—2)    Idy  (3.5) 

0     oyfa  l     '       Gy       J 
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If  we  let  z  =   (y-<H   )/o   ,   then  o  dz  =  dy.     When  y  =  0,   the  lower  limit 
of  z  is  given  by 


z  =  ^  = 


2 :  2 

°S  +  as 


(3.6) 


and  when  y  -  +  «9  the  upper  limit  of  z  ->  +  «. 
(3.5)  can  be  rewritten  as 


Therefore,  equation 


R--!-f 


-z 


e    dz 


2»   y$  -  ys 


(3.7) 


2  ,  2 

aS  +  cs 


Obviously,  the  random  variable  z  =  (y  -  u  )/o  is  the  standard  normal 
variable.  Hence,  equation  (3.7)  can  be  changed  to 


R  =  1  -  6 


\     yS  "  ys  ^ 


2  ,  2 

°S  +  as 


(3.8) 


Reliability  can  be  found  by  referring  z  = 
table. 


US   ■  ys 


2    .      2 
°S  +  as 


to  the  normal 


3.2     Reliability  computation  for  log-normally  distributed  strength  and 
stress  [10]. 
The  log-normal   density  function  is  given  by 


f  (p)   ■  — L-    exp  [-  -U  (*np  -  -,)2]  p  >  0 

M  pa/2TT  2o 


(3.9) 
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where  p  is  the  random  variable.     The  parameters  jj  and  a  are  the  expected 
value  and  the  standard  deviation,   respectively,  of  the  variable  £np 
which  is  normally  distributed. 
Let     x  =  £np.     Then  dx  =   (1/p  )dp.     From  equation   (3.9),  we  have 

f(x)   =  — —  exp  [-—■  (x  -  p)2]  .  co  <  x  <  -         (3.10) 

a/27T  2a 

V 

If    p    denote  the  median  of  p,   then  we  may  write 

v 

Pi  1  9 

0.5  =  /      — ■ — exp[-  —  (£np  -  ^)^]dp  (3.11) 

0        pa/2-n"  2a 

Using  the  transformation  x  =  £np,  we  rewrite  it  as 


£np         ,  ,  ? 

0.5  =  /  — —    exp  [-  -^  (x  -  y)^]dx  (3.12) 

-K         av2ir  2a 

v  V  II 

Yielding  u  =  £np       that  is,    p    =  e   . 

Returning  to  the  original    problem  in  which  S  and  s   are  log  normally 
distributed,     we  let  y  =  S/s,  which  means   £ny  =  £,nS  -  &ns.     £ny  is 
normally  distributed  since  both  £nS  and  £ns  are  normally  distributed. 

The  log  normal   density  function  is   positively  skewed  and  hence  the 
median  is  a  better  and  more  convenient  measure  of  the  central   tendency 
for  the  log  normal   distribution  than  the  mean.     The  antilogof  the  mean 
of  2,nS  is  the  median  of  f<.(0»   and  that  of  £ns  the  median  of  f  (•)>   that 
is, 
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y  _     uanS 
5  -  e 


or     ,mS  -  ml 


and 


v        ^ 
s  =  e 


ans 


or  u        =  ans 


where  b  and  s  are  the  medians  of  S  and  s   respectively,     By  analogy,  we 
add 


yany  =  W 


(3.13) 


since  the  variable  y  is  also  log  normally  distributed.  But, 


yany  =  yanS  "  yans  =  ^  "  *"* 


(3.14) 


Combining  the  two  equations,  (3.13)  and  (3.14),  we  get 


any7  =  anS  -  ans  =  an  y 

s 


and 


(3.15) 


2 


C  -       I  0  r-      +      G 

any      ;/    ^nS         ans 


(3.16) 


From  the  definition  of  reliability,  we  have 


R  -  P(|->  1)  -  P(y  >  1)  -  /  f(y)dy 
s  1     y 

Let     z  =   (any  -  n       )/on     .     Then  z  is  the  standard  normal   variable. 
v     J         any7     any 

To  find  the  new  limits  of  integration,  when  y  -  1, 


an  1    -  u 


z  = 


any  _  _    a 


nS  -  ans 


any 


r2      +     2 
aanS      °ans 


(3.17) 
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and  when  y  -»■  +  «,   z  ■>  +  ~.     Then,   the  reliability  can  be   rewritten  as 


R  "—  /     y       v  e"2  /2  d*  (3.18) 

/2tt  -£nS-£ns 


2         2 
£nS     £ns 


3.3  Reliability  computation  for  exponentially  distributed  strength 
and  stress  [10] 

The  pdf  of  an  exponentially  distributed  strength  S  is  given  by 

-ASS 

fs(S)   =  As  e  0  <_S  <  -  (3.19) 

and  that  for  stress   s  is 

-A    S 

f_(s)   =  a  e     s  0  <  s  <_  -  (3.20) 

Using  equation   (2.16),  we  have 


R  ■  /  fc(s)C/  fc(S)dS]ds 


00        -AS         °°        -A^S 

I   A  e     s     [/  AQe     ^     dS]ds 
0     s  s     b 


■     ,.        "Ac 

/  >.  e     s   [e       ]ds 
0     s 


.  A  .-Ws* 


0    s 
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A 
S 

h+xs 


(3.21) 


(3.22) 


If  the  mean  value  of  strength  is   denoted  by  p~,   then 


Uc  =  E(S)   =  /  f.(S)    •    SdS 
3  0 


00       -X-S 

/  >   e     b     •   SdS 
0     b 


=  1_ 

AS 

and  the  mean  value  of  stress  by  y  is 

M,  =  E(s)  =  y-  (3.23) 

s        As 

then 


yS 
R  =  — *-  (3.24) 

uS+Ms 

3.4  Reliability  computation  for  gamma  distributed  strength  and  stress 

[10] 

The  gamma  density  function  of  a  random  variable  x  is  given  by 

}nyn-l    -Ax 
f(x)   =    l     r(n) n   >  0,   A   >  0,     0  <  x  <  -  (3.25) 

where  A  is  called  the  scale  parameter  and  n  the  shape  parameter  and 
r(n)  is  the  gamma  function.  The  pdf  for  strength  is 


and  the  pdf  for  stress  is 


Using  equation   (2.21),  we  have,  as  before,   for  y  =  S-s 


R  =   f  f,(y+s)f.(s)ds 

0     3 


Let  v  =  s/y,   r  =  u/x,  and  dv  =   (l/y)ds,   then 


also  let 


I  ym+n-l  e-[(l+(l+r)v)x]ydy=         r(m+n) 
0 


(l+(l+r)v)       A 


then  we  have 


where  r  =  u/x.     If  let  w  =  rv/(l+(l+r)v) ,  we  have 
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fs(s)-7Ws     e 


X>0,     m>0,     0<S< 


(3.26) 


Vsj       r(n)  s      e 


u>0,     n>0,     0<s< 


(3.27) 


.m  n 

A    u 


r(m)r(n) 


0     0 


r  i    ,    xm-l   -x(y+s)  n-1    -us    .    .  Ir>  ooN 

J   iy+s)       e     w     ys       e         dsdy  (3.28) 


R  = 


,m  n 

x  u 


r(m)r(n) 


/y 


m+n-1   -(l+(l+r)v)yx 


dy   (  (l+v)m"1vn_1dv         (3.29) 
0 


m-1   n-1 


R  _  r'r(m+n)    ,  (1+v)"   'v 
r^r(n>  b  (H(l+r)v 


m+n 


dv 


(3.30) 
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dw  = 


(1+r)vr  J   dv= 


r u+r)vr  .     ..  vr 

iTWFTVT-   (1+(1+r)v)2j   dv"  (1  +  (l+r)v) 


dv 


R  '  r(m)r(n)  £  (1'w^       w       dw 

Four  special   cases  are  briefly  discussed  here, 

1 .     If  x  =  u  =  1 ,     then 


R  -  r(m+n)        ,     Q+v) 
r(m)r(n)  a       ,,+? 


0       (l+2v) 


m-1    n-l 

v         .     _  r(m+n  _ 

m+n    "    v  "r(iP)r(n")Q 


1/2 


(l-w)m"1wn'1    dw       (3.31) 


2.     If  m  =  n  =  1 ,  then  S  and  s   are  exponentially  distributed  with 


R  = 


r(2) 


r/l+r 


r(l)r(l) 


dw  = 


1+r      u+; 


(3.32) 


Equation  (3.32)   is   the  same  as  equation  (3.21). 

3.      If  m  f  1 ,  and  n  =  1,   then  the  strength  has  a  gamma  distribution  and 
the  stress  has  an  exponential       distribution.      In  this   case,  we  have 


.    _     r(m+l )      .r/|+r  m_i  1     m  \     m 


r(mMT)  }Q 


P"T  A 


(3.33) 


4.  If  m=l ,  and  n^l ,  then  the  strength  has  an  exponential  distribution 
and  the  stress  has  a  gamma  distribution.  In  this  case,  we  have 
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r/Cl+r) 

n  r(n+1 )    r/v       ;     n-1,        nr(n)   /   r  ,n  ]       ,   r  ,n  ,,  ... 

R=  r(1)r(n)  /  w      ^^WrW  <3'34) 

3-5     Reliability  computation  for  Weibull   distributed  strength  and 

stress.   [10] 

The  pdf  of  Weibull  distribution  has  been  mentioned  in  Section 
2-1.  The  pdf  for  the  strength  and  the  stress  are  given  by 

fS(S)   =  F"  ("r2)  exp   [-(— ^)     ],     SQ<S>~  (3.35) 

and 

Cc       S-S^    CS  S-So    "S 

f(s)=^(— -±)  exp  [-(— -2-)     ],     s0<s<-  (3.36) 

s  »s       8$  cs  o 

respectively.     The  probability  of  failure  given  in  equation  (2.20) 
wi  1 1    be 


R  =  P(S  <  s)   =  /     [1    -  Fs(S)]fs(S)dS 


s-s     s    &c    s-s    es-1  s-s    h 

=  /     exp[-(-^)     ]/(-^)  exp[-(— ^)     ]dS  (3.37) 

S  us  bS       JS  US 


where 


c  P 

F  (S)  =  1    -exp[-(-A)  S] 


es 


s-so  ps 

Let  q  =   (—2-) 
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then 


and 


6,     S-S      ~S    ' 
dq  =  -r  (-t-2")  dS  (3.38) 

es      s 


1/6. 
S  =   q  9S  +  S0  (3.39) 


then  we  have 


R  =  P(S  <  s)   =   /     e"q  exp[-[/  q       5  +   (-%-£-)]     ]dq  (3.40) 

0  cs  9s 


0      n  eQ     1/e-        S  -s„     6s 

R  =  P(S  >  s)  =  /     e_q  exp[-(/  q       5  +  (-^n     ]dq 

Ps  es 

3-6     Reliability  computation  for  extreme  value  distributions 

As  we  mentioned  in  chapter  2.     The  pdf  of  the  smallest  value 
distribution  is  given  by 

9n(x)   =  nf(x)[l-F(x)]n']  (3.41) 

and  the  cumulative  distribution  function  is  given  by 

Gn(x)   =  1    -   [1    -  F(x)]n 

where 

x 
F(x)  =  /       f(x)dx 

—00 

There  are  three  asymptotic  distributions   for  the  smallest  value. 
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They  have  the  following  form: 
Type   I 


x-x. 


Gn(x)  =  1   -  exp[-exp  (-g-5-)] 


<    X    <    +    <*>,       6>0 


(3.43) 


Type  II 


x-x. 


Gn(x)   ■  1 


exp 


[-(-  -r-^)     ]        —  <  x  <  xn      e  >  o,  e  >  o 


(3.44) 


Type  III 


x-x 


Gn(x)   =  1    -  exp[-i 


x     <_x<^,       e  >  0,  p  >  0 


(3.45) 
Each  type  of  these  distributions   arises  when  certain  conditions 
are  met.     Type  I  arises  when  the  underlying  density  function  tends  to 
zero  exponentially  as  x  -  +  °°.     On  the  other  hand,   if  the  range  of  the 
density  function  is   unbounded  from  below  and  if  for  some  e  >  0,   3  >  0 
we  have 


lim     (-xr   Fy(x)  =  6 


the   limiting  distribution  of  the  smallest  value  is  of  Type  II.     Type  III 

distribution  arises  when   the  following  conditions   are  met: 

1).     The  range  for  the  underlying  density  function   is  bounded  from  below 

(i  .e. ,  x  _>  xQ). 

2).  F  (x)  behaves  like  (e-x  )  ,  for  some  e  >  0,  e  >  0  and  as  x  ->■   ». 


The  pdf  of  the  largest  value  is  given  by 
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h  (x)   =  nf(x)[F(x)] 


n-1 


(3.46) 


and  the  cumulative  distribution   function  is  given  by 


Hn(x)   =   [F(x)]n 


(3.47) 


There  are  also  three  asymptotic  distributions  for  the  largest  value 
Type  I 


x-x. 


Hn(x)   =  exp[-  exp[-   (-r^)]] 


-oc<x<^,     e  >  0 


(3.48; 


Type  II 


Hn(x)   =  exp[-(^)     ] 


x  _>  x  ,  5  >  o,    e  >  o 


(3.49) 


Type  III 


X-  X       R 

Hn(x)  =  expH--^)  ] 


X  ^  X    ,    6   >   0,      5    >    0 


(3.50) 


Let  us   consider  the  case  where  the  strength  has  the  Weibull    distribution 
and  the  stress   has   the  Type  II   largest  extreme  value  distribution,  we 
have 


ss    S-S.  V1  S-S_  es 


fs(S)=^(^f)  eXp[-(^f}     ]         So-S< 


(3.51) 


and 


F  (s)   -  exp  [-   {-p-)     S] 

s 


s      .  s   <  oo,     e.    >  0,     3S   >  0     (3.52) 
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Substituting  these  two  equations  into  equation   (2.17)   gives 


S-s„  "£s     £c     S-S„  V  S-S     5S 

R  =  /     exp[-(-s-2-)       ]  -f-  (-r-^)  exp[-(— 2.)     ]  dS  (3.53) 


o 
Let 


S_  es  6S       eS  US 


S-S     6S 


q=   (~^) 


tiS 


then 

6Q  S-Sn  SS_1 

dq  =  t^  (-— -)            dS 

CS  CS 


and 

1/ 


S   '    eS  +  So  =  S 


Hence 


—  r 


oo  g        ]/ft  s    -s  s 

R  =    I     e"q  expL-Er^q     "S  +  -V5-]       ]   dq  (3.54) 

0  3s  "s 

Similarly,   the  reliability  in  other  case  can  be  derived  which  either 

strength  or  stress  has  extreme  value  distribution.     For  example: 

1).     Weibull   distributed  strength  and  Type  III   smallest  extreme  value 

distributed  stress: 


6   -6   £nq   -  S     PS 
R  =    I     e"q  exp[-(^-^- -)      ]  dq  (3.55) 

b  ds 

where 

q  =  exp[-(V^)] 
9s 
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2).     Type  I   smallest  extreme  value  strength  and  Weibull   distributed 
stress : 

00  e         1/3        6  q-  s 

R  =   /  e-q  exp[-exp  /  (q       s   -  -f-^)]  dq  (3.56) 

0  bS  cs 

where 

s-s     °s 

q   =   (— -0-) 
^s 

3).     Both  stress  and  strength  tend  to  be     Type  I   largest  extreme 

value  distribution: 

R  =    !  e"q  exp[-exp[^  Unq  -  -|— ^jjdq  (3.57) 

6  cs  aS 

where 

S-6 
q  =  exp[-   (— -*-)] 

CS 

There  still  have  many  combinations  of  different  strength  and 
stress  distributions,  which  we  did  not  derive  here.  But,  the  derivation 
is  quiet  similar  to  what  we  derived  here.  The  equations  we  derived 
before  almost  include  the  integration.  Those  integrations  cannot  in- 
tegrate analytically.  Hence,  a  numerical  integration  method  is  used, 
which  is  shown  in  Appendix  A. 
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CHAPTER  4.  TIME  DEPENDENT  STRESS-STRENGTH  MODELS  (S-S-T) 

In  Chapter  3,  the  stress-strength  interference  (S-S-I)  models  have 
been  presented.  However,  the  important  time-factor  is  not  considered. 
As  time  proceeds,  everything  changes.  A  physiological  or  psychological 
function  changes  as  a  human  ages.  A  radioactive  element  reduces  its 
weight  as  time  passes.  Concerning  society,  population  changes  when  time 
passes.  Similarly,  also  affects  the  stress  or  strength  distribution. 
The  time  factor  can  be  classified  into  the  following  two  categories:  [1-5] 
1).  Cyclic  Occurrences: 

The  random  variable,  such  as  stress  or  strength,  changes  according 
to  cyclic  pattern.  These  variations  may  be  due  to  natural  seasonal 
changes  or  day/ night  temperature  cycles  or  man  made  on-off,  up-down, 
cycles.  The  cyclic  changes  have  a  fixed  period  and  a  nonconstant  time 
between  cycles.  The  reliability  is  a  function  of  the  number  of  cycles, 
n,  rather  than  time.  Examples  are  the  operating  cycles  of  a  relay  or 
the  number  of  take  offs  and  landings  of  an  aircraft. 
2).  Random  Occurrences 

In  this  category,  the  times  between  variable  occurrences  are  random 
rather  than  known.  If  it  is  assumed  that  stresses  occur  randomly  and 
that  each  stress  occurrence  is  independent  of  each  other,  then  we  have 
a  Poisson  probability  law  for  stress  occurrence.  Random  occurrences  in- 
clude not  only  the  Poisson  probability  law  but  other  occurrence  laws 
as  well . 
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The  stress  and  the  strength  may  vary  during  a  long  operating  inter- 
val, in  relation  to  the  passage  of  time  or  to  the  number  and/or  severity 
of  previous  stresses.  We  will  use  the  following  classification  for  such 
time  variations. 

A).  Aging  -  Aging  describes  changes  with  time  in  the  parameters  of 
the  model.  Most  commonly  this  is  modeled  as  a  shift  of  the  mean  and/or 
variance,  e.g.,  a  linear  decrease  in  strength.  A  simple  example  is  the 
corrosion  in  a  liquid  cooling  system. 

B).  Cyclic  Damage  -  An  item  may  experience  a  change  in  its  strength 
as  the  device  undergoes  repeated  operating  cycles.  Thus,  the  strength 
density  is  a  function  of  the  number  of  cycles  n.  An  example  of  this 
phenomena  is  the  shortened  life  of  a  light  bulb  which  is  subject  to  many 
on-off  cycles  rather  than  allowed  to  burn  continuously. 

C).  Cumulative  Damage  -A  device  is  said  to  suffer  cumulative  damage 
when  its  decrease  in  strength  is  determined  by  the  magnitude  and  the 
number  of  previous  stresses.  An  example  would  be  air  leakage  from  a  space 
craft  due  to  meteorite  collisions  puncturing  the  skin.  Larger  meteorites 
making  larger  holes  creating  more  air  leakage  may  De  assumed. 

These  two  terms  cyclic  occurrences  and  cyclic  damage,  are  often 
confused  as  are  random  occurrence  and  aging.  The  similarity  between 
cyclic  occurrences  and  cyclic  damage  is  in  considering  cycles  n  rather 
than  time  period.  The  main  difference  between  them  is  that  the  cyclic 
damage  is  for  strength  only,  whereas  cyclic  occurrences  is  for  both 
strength  and  stress.  The  similarity  and  difference  between  random 
occurrences  and  aging  are  the  same  as  between  cyclic  occurrences  and 
cyclic  damage.  The  only  difference  between  former  and  latter  is  one  for 
time  period  and  another  for  cycles  n. 
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The  uncertainty  about  the  stress  and  the  strength  variables  may  be 
classified  in  the  following  three  categories:  1)  deterministic  stress 
or  strength,  2)  random-fixed  stress  or  strength,  3)  random-independent 
stress  or  strength. 
1).  Deterministic  stress  or  strength 

The  variable  is  either  constant  or  varies  in  some  known  predictable 
manner.  If  both  stress  and  strength  are  known,  the  failure  is  deter- 
ministic rather  than  probabilistic.  The  device  succeeds  if  the  strength 
is  greater  than  the  stress,  and  fails  if  the  strength  is  not  greater  th'a- 
the  stress.  Of  course  known  or  deterministic  implies  that  the  manu- 
facturing process  is  well  controlled  so  the  parameters  are  predictable 
or  a  simple  nondestructive  test  is  available  to  determine  the  parameters. 
This  control  rarely  exists  in  the  real  world,  but,  in  some  cases, 
acceptable  approximation  may  be  used. 
2).  Random-fixed  stress  or  strength 

We  are  interested  in  the  behavior  of  the  variable  with  respect  to 
time  or  cycles.  Strength  is  a  random  variable  at  any  particular  instant 
of  time.  It  is  assumed  that  enough  data  has  been  recorded  in  the  past 
to  determine  a  probability  density  function  for  strength  or  stress.  It 
is  also  assumed  that  any  test  to  determine  the  variable  precisely  is  too 
costly  or  destructive.  The  word  "fixed"  in  this  classification  refers 
to  the  behavior  of  the  random  variable  with  respect  to  time  and/or  cycles. 
This  means  that  the  random  variable  cnanges  in  time  in  a  known  manner. 
Let  f-  (S  )  be  the  pdf  of  the  strength  random  variable  S  at  the  initial 
time.  Then  the  strength  S(t)  at  any  instant  of  time  t  is  given  by 
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S(t)   =  SQp(t)  (4.1) 

Where  (|>(t)   is  a   known  or  given  function  that  describes   exactly  how 
strength  decreases  with  time.     For  example,   let 

<j>(t)   =   (1    -   0.0001    t)  (4.2) 

which  means  that  with  the  passage  of  eacn  unit  of  time,   strength  de- 
creases  byO.01%  of  the  initial   strength  S    .     Similarly,   let  S.    denote 
the  strength  for  the  kth  application  of  load,   and  let  us   assume  that    ' 
S^  is  a  function  of  S     and  the  occurrence  number  k.     Thus,  we  have 

S,  =  S   k-a  (4.3) 


'k      Jo 


or 


Sk  =  SQ  +  bk  (4.4) 

In  the  above  equations,  the  strength  decreases  with  the  number  of 

load  occurrences  for  a  >  0  and  b  <  0,  respectively.  If  increase  in 
strength  occurs  because  of  work  hardening,  we  have  a  <  0  and  b  >  0. 
3).  Random-independent  stress  or  strength 

In  this  category,  not  only  is  a  variable  value  sufficiently  unknown 
so  that  it  is  well  described  as  a  random  variable,  but  succesive  variable 
with  respective  to  time  and/or  cycles  are  so  unrelated  as  to  be  sta- 
tistically independent.  Observation  of  one  variable  value  gives  no  in- 
formation about  the  size  of  the  subsequent  value.  Successive  stresses 
are  generally  independent.  Strength  will  vary  randomly  and  will  be  in- 
dependently from  cycle  to  cycle  only  if  it  is  being  affected  by  other 
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environmental   factors,   such  as  temperature  and  vibrations,  which  are  in- 
dependent of  the  process. 

4.1     Reliability  computations   for  cyclic  occurrences 

The  reliability  of  stress-strength  interference  cases  without  con- 
sidering the  time  factor  are  presented  earlier  in  this   paper.     The  in- 
dividual  expressions  for  reliability  with  respect  to  time  are  developed 
here.     We  obtain  these  expressions  by  combining  the  three  levels  of  un- 
certainty classifications   for  both  stress  and  strength  and  two  types  of 
cycles  occurrences.     There  should  have  18  (=  3x3x2)   expressions  plus 
aging,  cyclic  damage  and  cumulative  damage.     Totally,   21    expressions 
should  be  developed.     But  not.  all   of  them  are  going  to  be  developed  here. 
Because  aging  and  cyclic  damage  are  special    cases   in  the  cyclic  occurrences 
and  random  occurrences  respectively,   the  reliability  of  random  occurrences 
will    not  be  developed  here  either.     Only  nine  cases  of  the  cyclic  occur- 
rences and  cumulative  damage  are  developed. 

In  the  cyclic  occurrences,  we   like  to  use  R  ,   the  reliability  after 
n  cycles   (the  probability  of  not  having  a  failure  on  any  one  of  the  n 
cycles),   rather  than  R(t),   the  reliability  at  time  t,   the  argument  t 
being  continuous,   because   the  former  can  be  converted  to  the  latter 
very  simply  when  cycle  times  are  deterministically  known.     For  example 

R(t)   =  Rn  tn  <  t  <  tn+1  n  =  1,   2,    ...  (4.5) 

where  t-    is  the   instant  in  time  at  which  the  ith  cycle  occurs. 
Case  1.     deterministic  stress  and  deterministic  strength  [10] 
Let  s.j   and  S«,   i   =1,   2,    ...,   n  denote  the  stress  and  strength, 


46 


respectively,  on  the  ith  cycle.     Then 

Rn  =  P[Er  E2,    ..,   En]  (4.6) 

where  R     =  reliability  after  n  cycles  of  occurrences 

E-   =  event  of  no  failure  occurs  on  the  ith  cycle 
Hence,  we  have 


if     s-    >  S-  for  some  i        1    <_  i   <  n 


Rn  = 


(l        if     s.j    <_  S-  for  all      i       1    <_  i   <_  n 


Case  2.     deterministic  stress  and  random-fixed  strength   [10] 

Let  s  ,  a  constant,  denote  the  stress  and  S,  denote  the  strength  on 

the  ith  cycle.     The, 

si  =  so  -  ai      *  =  1' 2-  ••• 

where  a^  >^  0  are  known  constants.  Further,  the  a-'s  are  assumed  to 

be  nondecreasing  in  time.  The  pdf  of  S  ,  f<.  (S  ),  is  assumed  to  be 

o 
known.  Then 


P^  =  P<snlSn) 


=  P(s^  <  S  -  a  ) 
v  o  —  o    n' 


=  P<S0i5o  +  V 


■  /  ,   fS  <So>dSo 
s  +a    o 
o  n 


But 


Rn  =  P[Er  E2,  ....  En] 


=  P[E] |E2,  E3,  ...,  En]  x  P[E2,  E3,  ...,  En] 


=  P[E1!E2,  E3,  ...,  En]  x  P[E2  |E3,  E4>  ....  Ep] 
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x  P[E 


>     <-  n  ' 


3'    ^4 


-    El 


=   P[E1|E2,   E3,    ...,    En]  x  P[E2;E3,    E4,    ...,   En] 


x  ...  xP^iyxPtg 


(4.7; 


All  but  the  last  term  (P[E  ])  in  the  R.H.S.  of  equation(4. 7)  are  1 's  be- 
cause of  tne  restrictions  on  the  a- 's  which  cause  the  strength  Sj  to  de- 
crease in  time.     Hence, 


Rn  ■  PLEn]  "  ■  fS  <VdS. 


(4.8) 


s  +a       "o 
o     n 


Let  Rn  n_i   be  the  conditional   reliability  for  the  nth  cycle  given  that  the 
device  has  survived  the  previous   (n-1)   cycles.     Then 


P[En,      ....     E  J 


Sn.n-1   =  P[EniET  E2'    '    "   En-1]  =  P[Er    ....   E^-,]  =  Rn_- 


(4.9) 


Case  3.     deterministic  stress  and  random-independent  strength    [10] 

Let  a  constant  s     denote  the  stress.     Let  f-   (S)   be  the  pdf  of  the 

bi 
random  variable  strength  Si   on  the  cycle  i.     Since  successive  values  of 
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S-   are  independent,  we  get 

Rn  =  P[Er   E2,    ...,    En]  =  P[E.,]  x  P[E2]  x   ...    xP[En]  (4.10) 


where 


PLE,]  =  P(sn  <  S.)  =   /     f-   (S)dS         i  =  1,  2 n 

srt       i 

o 

In  particular,   if  the  pdf  remains   unchanged  over  time,  that  is,   if 


fs   (S)   =  f     (S)   =    ...   =  f     (S)   =  fs(S)  (4.11) 

1  2  n 


then 

00 

Rn  =   (P^))"   =   [/     fs(S)dS]n 
so 

Case  4.     random-fixed  stress   and  deterministic  strength    [10] 

This  case  is  similar  to  case  2  but  the  roles  of  stress  and  strength  are 

reversed.     Let  s-   =  s     +  b-       i   =   1,   2,    ...   denote  the  stress   in  cycle  i, 

where  b-'s  are  known  nonnegative  constants,  nondecreasing  in  time.     Let 

a  constant  S     denote   the  strength,   and  the  pdf  of  s   ,   f     (R  ),   is 

o 
known.     Then 


Rn  ■  p[En] 


=  p[sn^Sn) 


=  P(s     +  b     <  S   ) 
v  o         n  —    o' 


=  P(soiSo-bn» 
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S   -b 

"    /  '    fc        (Sn)dSo  (4J2) 

0  so       u       ° 

Case   5.     random  fixed  stress  and  random-fixed  strength   [10] 
Let  the  stress  be  given  by 

s  •   =  s     +  b. 
i  o  1 

and  the  strength 

Si   =  S0  -  a,.  1   =  1,  2,    ... 

where     a- 's  and  b • ' s  have  the  same   restrictions  as   in  cases   2  and  4, 

respectively.     The  pdf's   f     (s   )   and  f~   (S   )   are  assumed  to  be  known. 

o  o 

It  is   required  that  the  stress   is  nondecreasing  and  the  strength   is  non> 

increasing.     Hence,   arguments   identical   to  those  in  case  2  would  yield 


Rn  ■  P^ 


=  p(sniV 


=  P(s     +  b     <  S     -  a   ) 
y   o         n  -    o         ny 


=  P(s     <  S     -  a     -  b   ) 
v   o  -    o         n         n; 


■=i\(So)^VVbnfs„(S°)dSodSo  U-13) 

0       o  0  o 


In  special  case,  when  a-  =  b-  =  0,  the  equation  (4.1;,)  can  be  recog- 
nized as  the  standard  expression  for  SSI  model. 
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Case  6.     random-fixed  stress   and  random-independent  strength    [10] 

Let  f£(s)   denote  the  pdf  of  a  random  variable  s.     The  successive  random 

strengths  S-j ,  S^,    . ..,   SR  are  independent,   and  identically  distributed 
with  pdf,   f<~(S),  by  assumption 
Thus 


Rn  -  P[E-| ,   Ep5    ... » 


The  E-  event  means  that  E-   (S-  >  s) 
Hence 


V 


Rn  =  p[(s1  >  s)  r\  (s2  >  s)r\  (s3  >  s)a...  n(sn  >  s)] 


=  P[min   (S-|  ,   S2,    .. . ,   Sn)   >  s] 


(4. 14) 


The  distribution  function  of  the  random  variable 

Smin  =  min(Sl'  V    ••"  V  (4.15) 

is  the  smallest  value  distribution  which  has  been  mentioned  in 
Section  2.1 . 


Gn(S)  =  !-[!-  FS(S)]' 


Now,  equation  (4.18)  can  be  rewritten  as 


R  =  PCS  ,n  >  s] 

n    L  mm    J 


4.15) 


Hence,  using  equation  (2.21),  we  have 


Rn  =  i  V!s)[l  -Fs(S)]"ds 


(4.17) 
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Case  7.  random-independent  stress  and  deterministic  strength  [10] 
This  case  is  an  exact  reversal  of  case  3.  By  exchange  between 
stress  s  and  strength  S,  we  have 

Rn  =  P^)  x  P(E2)x  ...  x  P(En)  (4.18) 

where 

P(Ei)  =  P(Si  <  SQ) 

So 
=  f   f,  (s)ds 
0    5i 

where     S     is  the  known   strength  and  f     (s)  represents   the  pdf  of  the 

random  variable  stress  s-   on  the  ith  cycle.     In  particular,    if  f     (s) 

i  s1 

=  fSi(s)  =   ...   ■  fc    (s)  =  fc(s), 
then 


s2  sn  s 


Rn  =  [H^)f  =  {/  °  fs(s)ds}n  (4.13) 

Case  8.     random-independent  stress  and  random-fixed  strength    [10] 
This   case  is  similar  to  case  6  and  we  have 

Rn  =  P[(s1   <  S)r\  (s2  <  S)A  ...f\  (sn  <  S)] 

=  P[max(s1,    ...,  sn)   <  S]  (4.20) 

Let  sm=v  =  Max  (si>s0,    ...,  sj.     Then  the  largest  value  distribution 
max  i      l  n 

°'fr  sn,n^  is  given  by 
max         3  J 


52 


where 


Fn(s)   -   [Fs(s)]' 


F   (s)  =  /     f  (s)ds 

5  q  5 


and  hence 


Rn=^     fs(S)[Fs(s)]n   dS 


(4.21) 


Case   9.     random-independent  stress  and  strength    [10] 


Let  f     (s)   and  f-   (S)   denote  the  pdf's  of  stress  s-   and  strength  S. 
i  i 


respectively,   in  cycle  i   =1,2, 
dependent,   then 


Since  s-'s   and  S- 's  are  in- 


Rn  =  P^n'   E„-t    ■••>   ^ 


(4.22) 


=  P(En)   x  P(En_1)   x   ...   x  P(E]) 


n    P(E. 
i=l        ] 


where  P(Ei)  =  P(si   <  Si ) 


/     f     (s)    I  f     (S)dS  ds 
0        i         s     5i 


In  particular,   if  f   (•)  and  fr(-)   do  not  change  with  time,   then 


R     =     n     P(E..)   =   (P(E])) 
11       i=l  ' 


(4.24) 
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where  P(Ei)   is  the  reliability  of  the  stress-strength  interference 
models  which  were  discussed  in  Chapter  3. 

4.2     Reliability  in  case  of  cumulative  damage 

As  we  mentioned  before,   cumulative  damage  of  a  device  is  when  its 
decrease  in  strength  is  determined  by  the  size  as  well   as  the  number  of 
previous  stresses.     In  the  cumulative  damage  models,   if  it  is  assumed  that 
the  cumulative  damage  weakens  the  strength  of  the  part  at  each  occurrence 
by  an  amount  proportional   to  the  applied  stress,   a  cumulative  damage 
law  based  on  the  sum  of  the  applied  stresses  is  given  by 

Sn  =  Si  -  X   CiSi  (4-25) 

wnere  c.  's  are  proportionality  constants. 

Let  us  formulate  the  situation  where  the  strength  changes  with 
time  or  the  number  of  load  occurrences.  Let  S-  denote  the  strength 
after  i  load  occurrences,  and  S-  be  modeled  as  a  deterministic  function 
of  the  initial  (random)  strength  S  and  i.  Thus,  we  have 

One  example  of  the  above  function  may  be 

S-j  =  SQ  o(i)  (4.26) 

where  c(i)  is  a  monotonically  decreasing  function  of  i ,  i  =  1 ,  2 

Let  us  assume  that  the  stress  s  ,  s-,,  s?,  ...,  s  are  independent  with 

distribution  function  F  ,  F  ,F  ,...,F  ,  respectively.  Then 

so   sl   s2      sn 
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Rn  =  P^(so  <  So)n>  (sl    *  V^  (s2  *  S2)A-../Msn   <  Sn)]  (4.27) 

Since   the  strength  S-    is  a  known  function  of  the  initial   strength  S   , 
equation   (4.27)  may  be   rewritten  as 

Rn  =  P[(so  K  So}^(sl   *  h(SV   1)}^  '••^(sn  *  h(s0'h))]  (4*28) 

If  the  initial   strength,  S   ,   is  between  S  and  S+dS,   the  probability  of 
success   (defined  by  event  E)   as  given  by  equation   (4.29)  may  be  ex- 
pressed by  the  following,  because  of  the  independence  of  the  stresses:. 


P[E!S   <  SQ   <  S  +  dS]  =  P(sQ  <  S)   P(s1    <  h(S,l)) 


P(sn  <  h(S,n)) 


(4.29) 


Hence,   the  reliability  R     is 


Rn  =   /     F     (S)    F     (h(S,l)) 
0         o  sl 


F.    (h(S,n))f  (S)dS 
bn  3 


(4.30) 


where  "MS)  is  the  pdf  of  the  initial  strength.  If  S^  =  SQ  y(i)  and 
all  the  stresses  are  independent  and  identically  distributed,  we  have 


Rn  =  /  fc(S)  {  n  F_(S<|>(i))}ds 
1   0       li=0 


(4.31) 


R  given  by  equation  (4.31)  may  be  approximated  as  follows.  Let  F, 

1 1  ■ 


=  1 


Fs  and  SQ  =  h(S     0) 


Then 


Rn  =   f  fc(S)   [  n  (1  -  F _(h(S,i)))]ds 

n       0     ^  i=0 


(4.32) 


Now,   the  following  inequality  exists, 
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n  (1  -  F.  (h(S,i)))>  1  -  I     F.  (h(S,i))  (4.33) 

i=0      si  i=0  si 

and  thus  the  lower  bound  on  R  is  given  by 

R  >  1  -  I  /  MS)  F  (n(S,i))dS  (4.34) 

n      i=l  0   b    si 

n 
If    F  (h(S,i))   1,  then,  the  lower  bound  given  by  equation  (4.34) 
i=0  si 

is  fairly  close  to  R  .  Now,  the  probability  of  failure  on  the  ith 

stress  appl ication  is 


Pf(i)  =  /  MS)  FQ  (h(S,i))dS  (4.35) 

t     0   a     si 

Hence,  R  can  be  expressed  as 

Rn  >  1  -  I     Pf(i)  -   exp  [-  I     Pf(i)]  (4.36) 

i=0  T  i=0  T 


The   relations  given  by  equation   (4.36)   are  close  approximations   to  R 


when 


I     Pf(1)«  1 

i=0     T 
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CHAPTER  5  OPTIMIZATION  IN  ENGINEERING  DESIGN 

In  the  previous  chapters,  the  emphasis  of  the  study  is  on  using  re- 
liability as  a  criterion  of  the  engineering  design.  However,  in 
real  life,  not  only  reliability,  but  also  other  criteria  are 
to  be  considered  in  the  design  work.  Generally,  the  following  three  cri- 
teria are  most  likely  to  be  considered: 

(1 )  the  rel iabil i ty 

(2)  the  cost 

(3)  the  weight 

Traditionally,  only  one  of  these  criteria  has  been  considered.  Tne 
optimal  solution  for  one  criterion  is  surely  not  an  optimal  solucion 
for  others.  In  most  of  the  cases,  the  solution  is  not  even 
satisfactory,  if  other  criteria  were  also  considered.  However, 
the  decision  maker  (DM)  always  wants  to  attain  more  than  one  objective 
or  goal  in  selecting  the  course  of  action.  One  of  the  major  reasons  for 
the  scarcity  of  multiple  objective  formulation  and  consideration  in 
literature  is  that  until  recently,  almost  all  the  solution  strategies 
developed  involved  a  single  objective  function. 

The  growing  tendency  to  incorporate  more  and  diverse  criteria  in 
a  system  design  and  the  resulting  difficulty  in  consolidating  such  cri- 
teria have  provided  an  impetus  to  the  development  of  multiple  criteria 
methodology  in  various  forms. 

Recently,  multiple  objective  analysis  has  been  applied  to  a  wide 
variety  of  problems  [6].  However,  no  paper  has  been  found  applying 
multiple  objective  decision  making  methods  (MODM)  in  the  field  of  re- 
liability analysis  and  design  policy  making. 
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In  this  chapter,  two  methods  of  MO DM  [6]  are  to  be  used  to  de- 
termine the  optimal  policy  of  a  particular  design  problem: 

(1 )  Nonlinear  goal  programming  [7] 

(2)  Sequential  mul tiobjecti ve  problem  solving  technique  [SEM0PS][14]. 
Section  5.1  states  the  problem  and  to  develop  the  models  for  three 
criteria.  Section  5.2  applys  the  nonlinear  goal  programming  to  the 
problem.  Section  5.3  applys  the  SEMOPS  method  to  the  problem. 
Section  5.4  discusses  of  results. 

5.1  Statement  of  problem  and  model  development 

A  particular  design  problem  is  presented  here.  The  model  for  three 
criteria  is  developed. 

5.1.1  Statement  of  problem 

A  critical  tensile  element  (Fig.  8)  is  to  be  designed.  The  load  P 
acting  on  the  element  is  a  random  variable,  which  is  known  and  tends  to  be 
a  normal  distribution  with  mean,  P,  and  standard  deviation,  a   .     The 
tensile  element  nas  a  circular  cross-section.  The  radius,  r,  is  a  random 
variable  because  of  the  manufacturing  tolerance.  The  ultimate  tensile 
strength  of  the  material  used  for  the  element  is  a  random  variable 
because  the  properties  of  the  material  vary.  Normally  distributed 
radius,  r,  and  ultimate  tensile  strength,  S,  are  assumed.  [10] 

The  problem  is  to  determine  the  following  variables: 

(1)  the  mean  ultimate  tensile  strength,  which  can  be  used  to  decide 
the  material 

(2)  the  standard  deviation  of  ultimate  tensile  strength 
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Figure  8     Circular  cross-section  shaft. 
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(3)  the  radius  of  the  element 

(4)  the  manufacturing  tolerance  fraction 

The  problem  has  the  following  tnree  criteria  which  are  to  be  con- 
sidered in  most  of  engineering  design: 

(1)  the  element  reliability 

(2)  the  total  cost  of  element 

(3)  the  total  weight  of  element 

In  determining  the  optimal  design  policy,  Kapur  [10]  selected  reli- 
ability as  the  criterion  and  used  tne  Lagrangian  multiplier  method. 

In  this  chapter,  two  methods  for  multiple  objective  decision  making 
are  used  for  obtaining  the  optimal  design  policy  based  on  three  criteria. 

5.1.2  Model  development 

The  problem  is  to  determine  (1)  mean  tensile  strength  (2)  standard 
deviation  of  tensile  strength  (3)  radius  of  element  (4)  the  manufacturing 
tolerance  fraction  based  on  the  following  criteria: 

(1)  the  element  reliability 

(2)  the  total  cost  of  component 

(3)  the  total  weight  of  component 

Notations: 

2 
A         =  the  mean  value  of  cross-section,   in 

b.       =  the  constants  of  the  ith  parameter 

C         =  the  total   cost,   dollars 

3 
D         =  the  density,   lb/in 

L         =  the  length  of  the  tensile  element,   in 
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P  =  the  mean  value  of  load,  lb 

R  =  the  reliability 

r  =  the  mean  value  of  radius,  in 

S  =  the  mean  value  of  ultimate  tensile  strength,  psi 

s  =  the  mean  value  of  tensile  stress,  psi 

W  =  the  total  weight,  lb 

Z  =  the  standard  normal  variable 

a    =  the  manufacturing  tolerance  fraction 

2 
o«    =  the  standard  deviation  of  cross-section,  in 

Cp    =  the  standard  deviation  of  load,  lb 

a  =  the  standard  deviation  of  radius,  in 

os    =  the  standard  deviation  of  ultimate  tensile  strength,  psi 

o    =  the  standard  deviation  of  tensile  stress,  psi 

The  three  criteria  which  are  mentioned  above  are  evaluated  as  follows 
(1)  The  element  reliability 

The  strength  and  stress  are  assumed  to  be  normally  distributed. 
From  equation  (3.3),  we  can  derive 

Z  »  -j    S   =   b   s  (5.1) 

Because  Z  value  is  closely  related  to  the  reliability  in  equation   (3.8), 
equation   (5.1)   is  used  to  represent  the  reliability  criterion  in  this 
problem  in  order  to    avoid  integration  during  the  computation 
process . 
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2 
The  tensile  stress,   s,   is  given  by  s  =  P/A  where  A  =  zr   .     Using 

-2 
the  Taylor's  series  approximation  it  can  be  derived  that  A  =  -nr     and 

a«=2irra  .       If  the  tolerance  on  the  radius  of  the  circuluar  cross- 
section  is  a  fraction  a  of  r,   then   3c     =  ar;  that  is,  o     =   (ct/3)r.     Now. 
using  equations   (2.36)   and   (2.37),  we  have 


s   =  P/A  =  P/^r2  (5.2) 

and 

2         2     1  2     P 

S  Y      A  A      ^ 

-         -2  2  2-4     2 

Substituting  the  value  of  A  =  tit     and  c.   =  4-  r     a  /9  into  equation   (5.3) 

yields 

2        o2  +    (4/9)a2P^ 

7T  r 

Substitution  of  equations  (5.2)  and  (5.4)  into  equation  (5.1),  we  have 


Z  •  .    S  7^     m  (5.5) 

/  2   ^  +  W9)g  1 

This  relates  the  several  variables  to  the  Z  values  directly  and  to  the 

expected  reliability  indirectly. 

(2)  The  total  cost  for  making  or  purchasing  the  element 

The  cost  of  the  element  depend  on  the  values  of  S,  c<-,  r,  and  a. 
The  cost  function  was  proposed  by  Kapur  [10] 
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b2      ^     -b6      b8 
C  =  t^S   +  b3cs   +  bgs  °  +  byz     °  (5.6) 


The  parameters  of  each  variable  were  estimated  by  using  regression  ana- 
lysis and  by  assuming  other  variables  are  fixed.  The  variables,  s  and 
c  ,  in  equation  (5.6)  are  slightly  different  from  the  variables,  r  and 
a,  of  this  problem,  but  they  are  still  correlated.  Hence,  equation 
(5.6)  can  be  rewritten  as 

b?         b.      br  bo 

C  =  t^S  c   +  b3(o  /S)  q   +  b5r  D  +  bya  Q  (5.7) 

The  second  terms   in  equation   (5.6)  and  equation   (5.7)  are  different, 
because  the  same  value  of  cv  will   not  cost  the  same  for  different  values 
of  S.     For  example,   the  cost  for  o~  =  4,000  is  different  when  S  =  20,000 
from  S  =  14,000.     Hence,   the  ratio  of  c<-  to  S  is  to  be  prefered  in  the 
second  term. 
Another  cost  function  is  proposed  as  follows: 

-b?  bA  bfi  ba 

C  =  t^S       •    b3(c$/S)   H   •   b5  r  D   •   bya  (5.8) 

The  difference  between  equations   (5.7)   and   (5.8)   is  the  assumption. 
Equation   (5.7)    assumes  that  each  variable  is  not  interrelated  with 
others  and  equation   (5.8)  assumes  that  they  are  interrelated.     For 
example,  with  a  higher  value  of  5,   usually  a  harder  material,   not  only 
does   it  cost  more  for  the  material,   but  also  it  costs  more  to  produce 
accurately. 

A  set  of  data   is   needed   to  estimate  the  parameters  of  cost  function 
by  using  a  nonlinear  models  program  (N0NLIN) (Appendix  B).     Those  data, 
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which  include  four  variables  and  related  cost,  can  be  collected  by  asking 
the  cost  analysis  department  in  a  factory.  The  set  of  emperical  data, 
which  is  shown  in  Table  2,  was  collected  from  a  mechanic  of  Dept. 
of  Physics,  K.S.U  and  a  machine  shop  in  Manhattan  area.  Using  NONLIN 
computer  package,  the  estimated  parameters  and  the  sum  of  squares  and 
iterations  for  equations  (5.7)  and  (5.8)  are  shown  in  Table  3.  From 
Table  3,  equation  (5.8)  is  selected  as  the  cost  function  here  according  to 
less  iterations  and  smaller  sum  of  square,  that  is 

C  =  0. 5726(S/1 000. 0)1-04.0.5726(as/S)"0'4116.0.5726(r)0'801. 0.5726. 

(10000. a)"0'407 

where  the  constants   in  the  parentheses  of  the  first  term  and  fourtn 
term  are  used  to  prevent  the  overflow  and  underflow  during  the  computing 
processing. 

(3)   The  total  weight  of  the  element 
The  weight  function  is   given  by 

W  =  Density   •    Volume  (5.10) 

-2 
The  volume  is  equal   to  n  r  L,  where  ~  and  L  are  known,  and  r  is  variable. 

Hence,  equation   (5.10)   can  be  rewritten  as 

W  =  Density   .    tt . r2 . L  (5.11) 

The  density  in  equation  (5.11)  need  to  be  estimated.  The  density 
is  mostly  related  to  mean  tensile  strength  S. 


Table  2     Data  for  estimating  cost  function 
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^v.       variables 

S 

°s 

r 

a 

cost 

sample^\^ 

ksi 

ksi 

in 

xO.001 

dollars 

no                ^\^ 

1 

85.0 

5.10 

1.0 

4.0 

20.0 

2 

100.0 

10.80 

1.0 

4.0 

40.0 

3 

85.0 

5.10 

2.0 

4.0 

30.0 

4 

180.0 

10.80 

3.0 

4.0 

106.0 

5 

40.0 

2.40 

1.0 

4.0 

15.0 

6 

40.0 

2.40 

3.0 

4.0 

30.0 

7 

85.0 

5.10 

3.0 

4.0 

50.0 

8 

180.0 

10.80 

3.0 

1.0 

145.0 

9 

180.0 

10.80 

3.0 

2.0 

180.0 

10 

180.0 

10.80 

10.0 

4.0 

270.0 

11 

65.0 

3.90 

1.0 

7.0 

18.0 

12 

65.0 

3.90 

2.0 

4.0 

25.0 

13 

65.0 

3.90 

2.0 

2.0 

29.0 

14 

85.0 

2.55 

1.0 

4.0 

27.0 

15 

85.0 

2.55 

2.0 

4.0 

35.0 

16 

200.0 

5.4 

1.0 

4.0 

60.0 

17 

180.0 

1.8 

1.0 

4.0 

90.0 

18 

180.0 

14.4 

1.0 

4.0 

30.0 

19 

180.0 

18.0 

1.0 

4.0 

25.0 

20 

85.0 

5.1 

1.0 

6.0 

16.0 

21 

85.0 

5.1 

3.0 

2.0 

68.0 

66 


M-   <u 

LO 

=3- 

O    4- 

o\ 

o 

(0 

m 

• 

E    3 

CO 

O 

3  cr 

CM 

o 

oo    en 

p"" 

J"^ 

w 

c 

o 

•t— 

■M 

to 

r*. 

*■—■*. 

rC 

>=?- 

CO 

s_ 

• 

a> 

:-0 

+-> 

— 

i — i 

"O 

c 

IC 

O 

r^* 

CM 

o 

- — "> 

CO 

CO 

"3- 

r>. 

JQ 

, — 

o 

.n 

LO 

o 

c 

Ol 

CM 

o 

CM 

r^ 

•r— 

i->- 

CO 

m 

+-> 

-O 

. 

TO 

**■ 

o 

3 

-r 

cu 

aj 

4-> 
<D 
■= 
rs 
J. 
ra 
Q. 

ai 


-Q 


<3" 
CM 

co 


o 


<sO 

lO 

CM 

CO 

r^. 

"M 

-D 

en 


in 

o 

<4- 


Z5 

CJ 


CO 


VX5 


CO  , — 


CM 


o 

CM 

^ 

O 
1 

o 

o 
<* 

i 

O 

CTi 

O 

O  r— 


O 


<T3 

rx 
o 

C_J 


o 

CO 
CM 

CM 


CM 

LT3 


CO 

CD 


Z3 

cr 


CO 


67 


The  density  can  be  written  as 

b2 
D  =  b-,S  (5.12) 

A  set  of  data,  in  Table  4,  which  include  S  and  D,  is  collected  from 
ASTM  standard  [1],  Using  the  NONLIN  computer  package,  the  parameters  can 
be  estimated,  that  are, 

D  =  0.0746(S/1000.0)0'3  (5.13) 

where  the  constant,  1000.0,  is  used  to  prevent  overflow  during  the  compu- 
ting processing. 
Then,  equation  (5.11)  can  be  rewritten  as 


W  =  0.0746  (S/1000.0)°'3Tif2  .  L.  (5.14) 


Constraints 


(1)  Limits  of  mean  strength 


bul  -  S  1  bLl 


(2)  Limits  of  standard  deviation  of  strength 


u2  —  S  —  L2 


(3)  Limits  of  mean  radius 

b  0  <  r  <  b.  0 
u3  -   —  L3 

(4)  Limits  of  tolerance  fraction 


6ci 


Table  4     Data  for  estimating  the  parameters  of  density 

Density  (Ib/cu  in) 
0.270 

0.285 

0.283 

0.310 

0.290 

0.335 

0.100 

0.163 

0.323 

0.317 

0.277 


S  (ksi) 

1 

65.0 

2 

85.0 

3 

82.0 

4 

120.0 

5 

95.0 

6 

180.0 

7 

43.0 

8 

55.0 

9 

165.0 

10 

100.0 

n 

70.0 
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bu4  i  a  ^  bL4 
(5)  Limits  of  ratio  of  standard  deviation  to  mean  strength. 

bu5  1  VS"  i  bL5 

These  constraints  are  obtained  according  to  the  ranges  of  each  variables 
in  the  data  set  in  Table  2  and  4. 

Numerical  example: 

Let  X-. ,  Xp,  x~,  x.  represent  the  variables  S,  o~,  r  a  respectively 
and  f-i(x_),  f«(x),  f3U)  represent  Z,  C,  W  respectively. 

The  following  numerical  values  are  assumed 

P  =  4000  lb,  g  =  100  lb 

P 

L  =  6.0  in, 
bul    =  200,000  psi,  bL1    =  40,000     psi 

bu2  =     1S>000  Psi'  bij)  =  ^'00°       psi 

bu3  =     10.00  in,  bL3  =   0.07     in 

by4  =     0.0070,  bL4  =   0.0009 

bu5=    o.io,  bL5  =  0.01 

5.2     Nonlinear  goal    programming  for  engineering  design  problem    [9] 

For  the  nonlinear  goal    programming  the  decision-maker   (DM)   has  to 
setup  tne  goal   and   rank  of  importance  for. each  objective  [9]. 
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For  each  equality  or  inequality,  the  negative  deviational  variable 
d"  and  positive  deviational  variable  d  are  introduced,  which  d",  d   >  0 
and  d"  •  dT  =  0. 

If    fjCx)  -  bj, 

< 

then 

> 

The  three  possibilities  (=)  may  be  acnieved  by  minimizing  a  linear 

< 

function  of  the  deviational   variables  as  shown  in  Table  5. 
5.2.1      Case  1 : 


Suppose  the  DM  states  the  following  objectives  in  order  of  priority 
Max     f^x) 
min     fo(x.) 
min     fo(xJ 
and  the  goals  are  setup  by  the  DM: 

f1(x)_^4.0       that   is     R  >  0.99996833 
f2(x)   <_  0.140 


f3(x)   <   0.00040 


Hence,  the  design  problem  can  be  formulated  as  follows 
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Table  5     Procedure  for  achieving  an  objective   [9] 

Objective  Procedure 

(a)  Equal   or  exceed  b-  minimize  dT 

(b)  Equal   or  be  less   than  b.  minimize  d- 

(c)  Equal   b.  minimize  6~.  +  d. 

b.    is  the  right  hand  side  of  inequality  i   or  equality   i 
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Find     x_  =    (x-i ,   Xp>   x^,   x4)  so  as  to  minimize 

a=[(d-+d-+d-+d-+d>d-+d;+d-+d-+d-0),    (d",).^).^)] 


subject  to 


G^      x]    -  40,000.0  +  d~   -  dj  =  0 

G?:     200,000.0  -  x1   +  d"   -  d*  =  0  (5.15) 


r    tvw»"uw-u     ai      u2     u2 

G3:      x2  -   2,000.0  +  d"   -  d^  =   0 
G4:      18,000.0  -  x2  +  d~   -  d*  =   0 


G5:     x3  -   0.07  +  d~  -  dt  =  0 


Gc:      10.0  -  x,  +  d"  -  d*  =  0 
5  3  6  6 


G7:     x4  -   0.0009  +  d~   -  d*  =   0 


Gg:      0.007  -  x4  +  dg  -  dg  =  0 


+ 


'a  . 


0.1    -  x2/x]   +  dg  -   dg  =   0 


G10:   X2/xl    "   °'01    +  d10  "   d!o  =   ° 

x,    -   1273.23/x^  , 

G,,:  '  J  +  d"     -  d,,   =  4.0 

11       / j         11  il 

/   2       1013.2  +  720434.1    x^ 

Vx2  +  4 

x3 

G12:    0.5726   (x^lOOO.O)1 ,04    .05726(x2/x] )"°-41 16    .O.d726(x3)0'801 
0. 5726(1000. 0-x4)"0,407  +  d~2  -  d*     =  0.14 
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G13:    0.0746(x1/1000.0)0'3   .18.85  X3  +  d~3  -  d|3  =  0.0004 

The  first  ten   goals  are  the  absolute  constraints  in  the  original   problem 
and  the  last  three  goals  are  the  three  objectives. 

This   problem  is  solved  by  using  the  nonlinear  goal   programming 
computer  package  of  Department  of  Industrial    Engineering,   K.S.U.    [7]. 

The  initial    starting  point:      x°  =   (42200,    3752.5,   0.20,   0.004). 

The  step  size:         ix°  =   (50.0,   10.0,   0.10,   0.001) 

After  32  iterations,   the  final   solution  is 

x*  =   (40650,   4055,   0.23,   0.007) 

a_*  =   (0.0,   0.0,    1.581,   0.22446) 

The  results  from  the  achievement  function  a_*  mean  that  the  first  priority, 
which  consists  of  zen  absolute  constraints,   has  be^n  satisfied  and  the 
second  priority,  which  is  Z  >_  4.0,   has  been  satisfied  too.     The  third 
and  fourth  priorities  have  not  been  satisfied,  which  have  1.681   and  0.22446 
far  from  their  goals  respectively.     The  final    values  for  three  criteria  can 
be  obtained  by  substituting  x*  into  three  criterion  functions,  we  have 

Z  value  =  4.0  that  is,   Reliability  =  0.999,968,33 

Total    cost  =  1 .821   dollars 

Total  weight  =  0.22486  lb 
The  optimal   design  policy  is  that  the  material   for  the  element  should  have 
the  ultimate  tensile  strength,   S  =  40650  psi,   and  the  standard  deviation, 
a~  =  4065  psi,   the  radius  of  the  element  is  0.23  in  with  manufacturing 


tolerance  fraction  0.7%.     The  results  give  the  reliability  of  0.999,968,33, 
the  cost  of  1.821    dollars,   and  the  weight  of  0.22486  lb. 

According  to  ASTM  standard  [1],   the  ultimate  tensile  strength  of 
Aluminum  Alloy  5052  is  similar  to  the  final    solution  of  S  =  40650  ps: 

5.2.2     Case  2: 


If  the  DM  feels  that  the  cost  of  one  component,   1.821    dollars,  which 
is  from  the  result  of  Case  1,   cannot  compete  witn  the  same  products  from 
other  companies,   then  the  DM  would  like  to  reduce  the  goal   for  cost  cri- 
terion to   1.8  dollars   and  put  in  a  priority  over  other  two  criteria.     The 
problem  becomes : 

Min  f2(x) 
Max  f1 (x) 
Min     f3(x) 

and  the  goals  are  set  up  by  the  DM: 

f2(x)  <  1.8 

f-,(x)  >_  4.0 

f3(x)  <  0.2 

Hence,  the  problem  can  be  formulated  as  the  same  as  (5.15) 

The  differences  are  that  the  orders  of  the  second  and  third  priorities 

of 
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a  =   [(d1'+d2+d3+d^d^+dg+d7+d>dg+d^0),(d12),(d^),(d13)] 
are  exchanged,   in  addition  to  the  goals  on  f«  and  f-,   that  is, 

G12:    0.5726    •   (xj/1  000. 0)1  -04  •  0.5726   •    (x2/x]  )"0'41 16  •  0.5726(x3)0'801 
0.5726-  (1000.0   .    x4r0'407  +  d~2  -  d|2  =   1.8 

G,3:   0.0746-(x1/1000.0)0,3   .    18.35   .   x^  +  d~^  -  d*     =  0.2 

The  nonlinear  goal  programming  problem  was  solved  by: 

The  initial  starting  point:  x°  =  (42200,  3752.5,  0.20,  0.004) 

The  step  size:     ax°  =  (50.0,  10.0,  0.10,  0.001) 

After  43  iterations,  the  final  solution  is 

x*  =  (426^5.5,  3785.0,  0.2,  0.007) 

a*  =  (0.0,  0.0,  1.211,  0.0) 

The  results  from  the  achievement  function  a*  mean  that  the  goals  of  the 
first,   second  and  fourth  priorities  have  been  satisfied.     The  third 
priority,  which   is  Z  >_  4.0,   has  not  been  satisfied,  which  has  1.211   far 
from  the  goal.     The   final   values  for  three  criterion  functions  can  be  ob- 
tained by  substituting  x*  into  those  criterion  functions,  we  have 

Z  value  =  2.789,   that   is,   Reliability  =  0.99730 
Total   cost  =  1 .8  dollars 
Total  weight  =  0.17  lb 
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The  optimal   design  policy  is  that  the  material   for  the  element 
should  have  the  ultimate  tensile  strength,  S  =  42625.5  psi ,  the 
standard  deviation,   cs  =  3785  psi,   and  the  radius  of  element  is  0.2  in 
with  manufacturing  tolerance  fraction  0.7%.     The  element  reliability  is 
equal   to  0.99730,   the  cost  is  1.8  dollars,   and  the  weight  is  0.17  lb. 

According  to  ASTM  standard  [1],   the  ultimate  tensile  strength  of 
Aluminum    5086  is  similar  to  the  optimal    solution  of  §  =  42625.5  psi. 

5.2.3  Discussion: 

The  comparison  of  the  results  of  Case  1  and  Case  2  is  given  in 
Table  6.  The  results  are  different  between  two  cases,  because  the 
ranking  of  importance  and  goal  of  each  criterion  which  have  been  selected 
by  the  DM  are  different.  It  can  be  concluded  that  the  nonlinear  goal 
programming  problem  is  very  sensitive  to  the  ranking  of  importance  and 
goal  of  each  criterion  function.  From  Table  6,  the  values  of  X3,  that 
is  the  radius  of  element,  is  reduced,  the  value  of  three  criterion 
functions  are  also  reduced.  It  can  also  be  concluded  that  the  radius 
is  the  most  sensitive  variable  in  this  two  cases. 

5.3  Sequential  Mul tiobjecti ve  Problem  Solving  technique  (SEM0PS)  for 
engineering  design  problem 

5.3.1  The  method 

A  Sequential  Mul tiobjecti ve  Problem  Solving  technique  (SEM0PS),  pro- 
posed by  Monarchi,  Ki  si  el ,  and  Duckstein  [14],  allows  the  DM  to  trade 
off  one  objective  vs  another  in  an  interactive  manner.  Semops  cyclically 
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uses  a  surrogate  objective  function  based  on  goals  and  the  DM's  aspir- 
ations about  achieving  the  objectives.  The  goal  levels  are  conditions 
imposed  on  the  DM  by  external  forces  and  the  aspiration  levels  are  the 
attainment  levels  of  the  objectives  that  the  DM  personally  desires  to 
achieve.  (Goals  do  not  change,  but  aspiration  levels  change  as  each 
iterative  cycle  goes.) 

Let  A  =  (A-, ,  . . . ,  A,  )  be  the  DM's  aspirations  levels,  and 
f(_x)  =  (f-,(x_),  ...,  f^U))  be   the  multiple  objective  functions.  Five  types 
of  objectives  and  the  corresponding  dimensionless  indicator  of  attain-' 
ment,  d_,  are: 

(1 )  at  most: 

f.  (x)  <  Ai ;    di  =  f i  (x)  /  A. 

(2)  at  least: 


fi  (x)  >  Ai ;    d.  =  A.  /  f i  (x) 


(3)  equals: 


^(x) 


A,s 


if  i  (21) 


Ai 


(4)  within  an  interval  : 


AiLifiWi'i«;   di 


A, 


U If  HiL      |    fM]) 


iL  +  A.^tfTTzr    Aiyi 


(5)     outside  an  interval 


f,-W  1  AiL  or  f.j(x)  >  Ai(J;  d. 


fAiL  +  AiU'if 


AiU       HAiL 


1  1 


W  +  ^ 


1    — ' 


u 
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Types  (1),  (2)  and  (4)  are  the  most  common.  In  each  instance,  values 
of  d-  £  1  imply  that  the  objective  is  satisfied.  Except  for  the  first 
type  (at  most),  the  d-  are  all  nonlinear  functions  of  an  objective 
function  that  may  itself  be  nonlinear. 

The  algorithm  generates  information  under  the  guidance  of  the  DM 
so  that  he/she  can  make  a  decision.  Information  concerning  the  inter- 
relationships between  objectives  is  in  terms  of  how  achievement  or  non- 
achievement  of  one  objective  affects  the  aspiration  levels  of  other 
objectives. 

The  cyclical  optimization  of  a  surrogate  objective  function  s, 
(5.16),  is  the  mechanism  by  which  information  is  generated  for  the  DM. 
The  word  surrogate  is  used  in  recognition  of  the  fact  that  the  true 
preference  function  of  the  individual  is  unknown.  Let  T'  be  the  sub- 
set of  the  set  of  T  objectives  as  those  objectives  making  up  s  at  a 
given  iteration  of  the  decision-making  process.  Thus, 

s  =     d.  (5.16) 

teT' 

is  defined  as  the  surrogate  objective  function.     The  value  of  each  d     in 
s  reflects  whether  the  t-th  objective  has  been  satisfied;   unsatisfied 
objectives   have  values   >  1. 

Operationally,  Semops   is   a  three-step  algorithm  involving  setup,   iter- 
ation,  and  termination.     Setup  involves   transforming  the  original   problem 
into  a  principal   surrogate  objective  function  problem  and  a  set  of 
auxiliary  problems   involving  surrogate  objective  functions.     The  "Iteration 
step  is  the  interactie  segment  of  the  algorithm  and  involves  a  cycling 
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between  an  optimization  phase  (by  analyst)   and  an  evaluation  phase   (by  the 
DM)   until   a  satisfaction  is   reached,  which  terminates  the  algorithm. 

The  first  iteration,  i=l,  solves   the  principal    problem,  and  a  set  of 

T  auxiliary   problems  formed  as  follows,  where  the  aspiration  level   of  each 
objective  is  given  as  the  goal   of  each  objective,   i.e.,  A-   =  b., 
i   =  1,2,    ...,  T. 

The  principal    problem: 


mm 


T 
t=l      ] 


(5.17) 


The  set  of  auxiliary  problems,   A  =  1,2,    ...,  T: 


min         s,„  =     T     d. 

U  >,  t 


tf 


s.t.       x  e    X 


V*>>AX 


(5.18  ) 


Solving  (5.17)  and  (5.18)  forms  the  optimization  phase.  The  re- 
sulting policy  vector  and  objectives  for  the  principal  problem  and  the 
set  of  auxiliary  problems  are  presented  to,  and  are  used  in  the  evaluation 
phase  by  the  DM.  The  impact  of  an  action  on  the  attainment  of  the  other 
objective  is  assessed,  and  a  new  aspiration  level  for  an  objective  is  set. 

In  general,  the  ith  iteration  solves  the  following  principal  problem, 
and  a  set  of  auxiliary  problems. 
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The  principal  problem 


min  s-  =     d. 
1   teT' 


s.t.  x  e  X 


(5.19) 


fj(x)  >  Ajs    j  e  (T-T1) 


The  set  of  auxiliary  problems,   HeT1    (the  number  of  T1   =  T  -   i+1) 


mm     s. 


=   y    d 

til 


s.t.   x  c  X 


fj(x)  >  Aj         for     Vj5       j    e   (T-T1) 


f   (x)   >  A 

X.    —     —      : 


ror  one  l,  I .  e  T 


(5.20) 


The  optimization  pnase  solves  (5.19)   and  (5.20).     The  resulting  solutions 
are  used   in  the  evaluation  phase  and  a  guidance  is   given  by  the  DM  for 
the  next  iteration  cycle. 


5.3.2     Tne  design  problem: 

Tnis  study  of  finding  the  optimal   design  policy  involves  three 
goals   and  four  decision  variables.     Their  various  formulations  are 
expressed  below: 
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Goals  and  criterion  functions 


x,    -  1273. 23/x^ 

fl  W  =  '  -  1  Al 

HJ 1013.2  +  720434.1    •    x< 

x2+  4 

X3 


f2(x)   =   0.  5726(x1/1000. 0)1'04   .  0.  5726(x2/x1  )"0>41  ]  6   .0.  5726(x3)0'801 

0.5726(1000.0  -  x4)"0,407  <  A2 


f3(x)  =  0. 0746(x1/1000. 0)0,3  •    18. SSx^  <  A3 


Goal    levels    (initial   aspiration   levels) 

GL-j    =  4.0,  GL2  =   1.7  dollars,   GL3  =   0.2  lb 

wnere     GU    =  4.0  is  equivalent  to,    Reliability  =   0.99996833. 
Relevant  range  of  f,-(x_) 

r(f1(x))  =  [-40.0,  lo.o],  r(f2(x))  =  [0.14,  10.0] 

r(f3(x))  =  [0.0005,  10.0] 


Constraint 
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200,000  <_x]    <  40,000 
18,000   <x2i     2,000 


10.0  <  x3  <_     0.07 


0.007  <x.<  0.0009 


0.1   <  x2/x]   <  0.01    j 


£(x)    <_  0 


Set  up: 

In  order  to  avoiding  the  negative  value  of  d-,   the  criterion 
functions  and  the  aspiration  levels  are  needed  to  be  normalized  as  follow: 


Vii)   "  fj(x)L 
Yi   =  fiWu-  fi(x)L 


AN.   = 


Ai    -fi«L 


1       fi^)u  -  fi(x)| 


where  "T -  (>S_) i  anc^  ^i^— \i  are  tne  ^ower  an^  upper  relevant  range  of  f  •  (xj 
Hence,  the  three  criterion  function  are  normalized: 


Yi  = 


f^x)  +  40.0 


50.0 


Y2  = 


f2(x)  -  0.14 
9786 


f3(x)   -   0.0005 
Y3  =  9.9995 


The  initial   aspiration  levels   are  assumed  to  be  equal   to  the  goal   levels, 
A.j    =  GL.j,    i   =  1,2,3  and   so  the  values  of  AN.   are: 
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AN1=M=0-88'     AN2  =  of=o-1582 


AN3  =  §7§lfi-=  0.01995 


The  three  indicators  of  attainment  are 


d1  =  AN1/Y1 


d2  =  Y2/AN2 


d3  =  Y3/AN3 


The  first  cycle: 

The  principal  problem  to  be  solved  on  the  first  cycle  is 


min   s-j  =  d-,  +  d2  +  cL 


s.t.  £(x)  <  0 

Also  as   part  of  the  first  cycle,  we  construct  three  auxiliary  problems 
which  attempt  to  satisfy  each  of  the  goals   in  turn.     If  f -,  (xj  ^4.0  is 
entered  as  a  constraint,   then  d-,    is  deleted  from  the  surrogate  objective 
function  givings 


min     S-,    -I    =  cL  +  <jL 


s.t.  £(x)   <_  0 


f^x)   >  4.0 


Similarly,  the  second,  and  third  auxiliary  problems  are: 
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min    S,  2  =  d-i  +  d3 


s.t    gjxj  ^  0 


f2(x)  il.7 


min 


S1.3  =  dl  +  di 


S.t     £(x)  1  0 


f3(x)  i  0.2 


Each  problem  is  solved  by  using  Sequential  Unconstrainted  Mini- 
mization Technique  (SUMT)  computer  package  of  Hwang  et  al .  [8].  The  re- 
sults for  the  principal  problem  and  three  auxiliary  problems  are  shown 
in  Table  7.  From  Table  7,none  of  the  results  is  satisfied.  Hence,  the 
procedure  goes  into  the  evaluation  phase.  On  examination  of  those  numbers, 
it  is  apparent  that  the  change  of  goal  2,  which  is  cost,  is  relatively 
independent  of  attainment  or  nonattainment  of  the  other  goals.  Also,  the 
DM  feels  that  he/she  would  like  to  raise  the  goal  2  level  (cost)  in  order 
to  get  higher  value  of  goal  1  level.  So,  it  seems  reasonable  to  choose 
an  aspiration  level  for  goal  2  and  enter  it  as  a  constraint. 

Using  Table  7,  the  DM  can  assess  the  impact  of  such  an  action  on 
the  attainment  of  the  other  goals.  This  assessment  is  made  as  follows. 
The  results  of  the  principal  problem  and  the  second  auxiliary  problem  are 
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Goal  Principal    Problem  Auxiliary  Problem 

1  -7.545  -3.794 

2  1.408  1.670 

3  0.0727  0.0929 

where  goal  1  is  pure  number,  goal  3  is  the  total  weight  in  lbs. 

The  new  aspiration  level  for  goal  2  on  the  next  cycle  is  Al.  Let 
q  be  the  proportion  (A«  -  1 .670)/(l .408-1 .670) .  Then  the  approximate 
effects  of  adding  these  constraints  are 

Goal  Approximate  Effects 

1  q  •  (-7.545  -  (-3.794)  +  (-3.794) 

3  q  •  (0.0727  -  0.0929)  +  (0.0929) 

The  DM  can  try  out  several    values  of  A?  until   he  believes  that  he  can 
accept  the  estimated  effect  on  the  other  goals.     This   procedure  estimates 
the  results  for  the  principal   problem  of  the  next  cycle.     It  assumes 
that  df-i  (x)/df?(x)   and  df_(x_)/df?(x_)   are  constants,   that  means   f-.(x)  and 
f2(x_)  or  f3(x_)  and  fo(x.)   are  linearly  related.     The  DM  sets  a  new  aspir- 
ation level   of  2.1   for  goal    2  and  estimates   the  results  as 
Goal  Results 

1  2.361 

3  0.1261 

where  goal    1    is  pure  number,   goal    3  is  the  total  weight  in  lbs.     The  DM 

is  satisfied  with  these  new  aspiration  levels  and 

to  2.1   for  goal   2  is  not  greatly  different  from  original   goal   of  1.7 


The  second  cycle 

The  new  aspiration  level  A  =  (4.0,  2.1,  0.2)  and  goal  2  is  entered 
into  constraint.  Proceeding  to  the  second  cycle,  only  two  auxiliary 
problems  are  necessary  solve  because  one  goal  has  been  added  to  the  con- 
straint set  for  this  principal  program. 
The  principal  problem: 

min    s2  =  d-,  +  cL 
s .  t .   £(  x )  i  0 

f2(x)  <  2.1 

The  auxiliary  problem  1: 


min         s~  -J   =  cL 


s.t.       £(x)  <  0 

f^x)   >  4.0 
f2(x)   <  2.1 
The  auxiliary  problem  2: 


min        s?  3  =  d, 


s.t.       f2(x)  <  2.1 

f3(x)   <_  0.2 

■-(■-s)'<  n 
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Each  problem  is  solved  by  SUMT  computer  package  of  Hwang  et  al  [8].  The 
results  for  the  principal  problem  and  two  auxiliary  problems  are  shown  in 
Table  8.  From  Table  8,  the  auxiliary  problem  1  does  not  have  a  solution, 
because  the  feasible  region  is  empty.  The  rest  of  results  show 
that  none  of  results  is  to  be  satisfied,  but  goals  2  and  3  have  been 
satisfied.  So,  the  DM  decides  to  decrease  the  aspiration  level  of  goal  1  ; 
that  means  the  DM  reduces  the  Z  values,  reducing  the  element 
reliability.  Finally,  the  DM  decides  to  set  a  new  aspiration  level  of  2.5 
for  goal  1  and  add  it  to  be  constraint.  The  DM  argues  that  a  new  aspir- 
ation 2.5  for  goal  1,  which  corresponds  to  the  reliability  of  0.99379, 
can  be  accepted. 
The  third  cycle 

The  new  aspiration  level  A  =  (2.5,  2.1,  0.2)  and  both  goals  1  and  2 
are  entered  as  constraints.  The  principal  problem  is: 

min    s-  =  d_ 

s.t    £(x)  <  0 

f](x)  >  2.5 

f2(x)  <  2.1 

and  there  is  no  auxiliary  problem.  The  problem  is  solved  by  using 
SUMT(LAI)  computer  package  [8].  The  results  are  tabulated  in  Table  9. 
From  the  Table,  the  DM  decides  that  he  is  content  with  the  results  and 

terminates  the  procedure. 
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5.3.3  Discussion 

The  results  in  Table  9,  which  are  satisfied  by  the  DM,  are 

x1    =  45772,     x2  =   3102.5,   x3  =  0.184,     x4  =  0.0064 

and     f^x)   =  2.515       f2(x)   =  2.09,     f3(x)  =  0.15 

The  results  mean  that  the  material  should  have  the  mean  ultimate  tensile 
strength,  S  =  45772  psi ,  and  standard  deviation,  cs  =  3102.5  psi,  and  the 
radius  of  the  element  is  0.184  in  with  manufacturing  tolerance  fraction 
0.64;;.  The  reliability  is  equal  to  0.99379,  the  cost  is  2.09  dollars,  and 
the  weight  is  0.15  lb.  Alluminum  alloy  5456  is  selected  for  the  material. 

5.4  Discussion 

The  two  methods  of  M0DM  are  used  to  solve  the  design  problem  in  this 
chapter.  Two  cases  of  nonlinear  goal  programming  are  studied  in  the 
Section  5.2.  The  results  of  two  methods  are  tabulated  in  Table  6.   The 
different  results  between  the  two  cases  in  the  nonlinear  goal  programming 
are  discussed  in  Section  5.2.3.  It  is  noted  that  the  results  are  in- 
fluenced by  the  DM's  ranking  and  goal  of  each  criterion.  The  solution 
is  very   sensitive  to  the  goal  vector  set  for  the  criterion  and  ordinal 
ranking  given  by  the  DM. 

The  optimization  technique  used  in  SEM0PS  does  not  solve  problems; 
it  generates  information  so  that  the  DM  can  look  at  the  inconsistent  con- 
straint set  and  selected  an  acceptable  alternative.  The  key  concept  of 
SEM0PS  is  its  interactive  nature.  This  prevents  the  problem  of  specifying 
the  DM's  preference  structure  by  allowing  him  to  keep  within  himself  the 
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transformation  that  he  makes  to  convert  numbers  into  value  judgements. 
The  DM  can  develop  a  ranking  of  goals  as  he  receives  information  con- 
cerning the  feasible  alternatives.  He  may  also  revise  his  preferences 
during  the  course  of  the  iterations.  In  the  design  problem,  it  shows 
that  SEMOPS  prompts  the  modification  of  conflicting  aspiration  levels 
so  that  an  acceptable  solution  can  be  determined.  SEMOPS  accomplishes 
this  by  revealing  to  the  DM  the  extent  to  which  his  aspirations  will  have 
to  be  modified  to  achieve  a  feasible  alternative. 

Two  methods  have  been  discussed  above.  The  preference  of  the 
methods  is  still  based  on  the  assumption  and  concept  of  each  method. 
If  the  DM  has  some  preference  and  goal  for  each  criterion,  then  the 
goal  programming  is  selected  to  use.  Otherwise,  the  SEMPOS  is  preferred. 
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CHAPTER  6  CONCLUSION  AND  DISCUSSION 

This  study  discusses  the  relationship  between  reliability  and  safety. 
It  also  discusses  several  models  for  computing  the  component  reliability. 
The  main  point  of  this  study  is  to  develop  the  design  model  which  is  based 
on  three  criteria: 

(1)  the  component  reliability 

(2)  the  total  cost  for  making  or  purchasing  the  component 

(3)  the  total  weight  of  the  component 

To  determine  the  optimal  design  policy  for  finding  (1)  the  mean  ulti- 
mate tensile  strength  (2)  the  standard  deviation  of  strength  (3)  radius 
of  the  component  (4)  the  manufacturing  tolerance  fraction,  which  satisfy 
the  three  criteria,  two  methods  of  multiple  objective  decision  making 
methods  are  used.  They  are  nonlinear  goal  programming,  and  the  Sequential 
Mul ti objective  Problem  Solving  technique  (SEMOPS).  The  advantages  of 
goal  programming  are  that  the  DM  does  not  need  to  give  the  numerical 
weights  for  the  objectives,  he  only  needs  to  give  an  ordinal  ranking  of 
them.  The  disadvantages  is  that  goal  programming  is  very  sensitive  to 
the  goal  vector  set  for  the  objectives  and  ordinal  ranking  given  by  the 
DM.  Goal  programming  has  been  widely  used  in  many  MODM  problems.  Se- 
quential Multiobjective  Problem  Solving  techniques  (SEMOPS)  belongs 
to  an  interactive  method.  The  advantages  of  this  method  are: 
(1 )  it  can  be  used  to  solve  nonlinear  problems,  and  (2)  the  DM 
can  reevaluate  his  desirable  achievement  levels  for  the  objectives. 
The  disadvantages  are:   (1)  since  the  optimization  phase  uses  an 
algorithm  based  on  linear  programming;  the  amount  of  linearization  work 


95 


is  quite  large;  (2)  it  is  handicapped  by  the  requirement  of  differenti- 
ability of  objective  functions  and  constraint  functions.  The  disadvan- 
tages may  be  eliminated  if  a  proper  nonlinear  programming  method  is  used 
such  as  SUMT  method. 

The  model  presented  in  this  article  only  restricted  to  three  ob- 
jectives; reliability  function,  cost  function,  and  weight  function.  But 
in  the  real  world's  design,  the  designer  needs  to  consider  more  criteria. 
Hence,  more  objectives  will  be  involved  such  as  the  reliability  for  shear 
stress  and  strength  interference  model  or  for  compressive  stress  and 
strength  interference  model.  How  to  estimate  the  cost  according  to  these 
parameters  will  need  further  studies. 

In  the  model,  one  of  the  objectives  concerns  reliability.  We  use 
Z,  the  standard  normal  variable,  to  represent  the  reliability  because 
the  integration  in  the  reliability  model  only  can  be  solved  by  numerical 
integration,  which  needs  too  much  computing  time. 

In  this  thesis  the  reliability  model  presented  belongs  to  an  SSI 
model.  The  SST  models  are  not  considered,  due  to  the  following  reasons: 
a)  Some  of  SST  models  rarely  happen,  but  they  are  almost  the  same  as 
SSI  models,  if  we  can  solve  SSI  models,  there  is  no  problem  to  do 
SST  models,  b)  although  some  of  the  SST  models  are  more  realistic  than 
some  SSI  models,  their  solutions  are  too  costly  in  computer  time,  and 
frequently  inaccurate. 
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APPENDIX  A 
The  Numerical  Integration  Program 
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COMPILER  OPTIONS  -  NAME  =   M  A  IN,  OPT*C2  ,  L  I  S'ECNT  =  60  .  S  IZE=COOO*  . 

S3URCE,E8CDIC,N0LIST  ,  NC2EC*  ,  IOAC.  "AP  ,  NOED  I  T  ,  ID.NOXREF 

Q  a  a  *  a  a  a  a  a  a  •  a  a  a  a«  «  *  a  »  c  r.  a  <t  aca  ae 

c  •  * 

C      *   NUMERICAL   INTEGRATION   « 

C      •  * 

£         a  a  a o a  a  a » a  a  *  a  a • c a  •  3  a  a  a  a  a a  »  a • a 

IMPLICIT     SEAL*8 t A-H.O-ZJ 
DIMENSION    XI 50) 
RE40I5, I) (XI [), 1-1,4] 

1  FORMAT (4F  10.21 
Nl-100 

Z*-l.*( XI 11-1273. 23 /(X(3>*«2)) /OSORT (  (X(2)*«2J*  ( (1013.2*720434 .l*( 
1X(4)**?  )  I /(X (3)**A  )  )  ) 

IF(Z)40,50,60 
AC    RX1=0EXP(-(Z*«2 )/2) 

RXQ«l.O 

H1»(0.-Z)/N1 

M=Nl/2 

SUH=RX0+RX1 

C015    J»1,M 

VVRR=0.0-M2*J-1.0)*Hl 

RX*OEXp(-(VVRR**21/2> 

SUM*SUM+4*RX 
IS  CONTINUE 

IF<H-l)26,26,13 
12  Ml=M-l 

CG25  J»1,M1 

WRR»0.C*>2*J»H1 

RX=DEXP(-( VVRR«*21/2J 

SUM=SUM+2  *RX 
25  CONTINUE 
2t    TN«Hl/3.* SUM/2. 5C66 

TN»TN*0.5 

GO  TO  IOC 
6C  RXl=OEXPl-(Z**2)/2I 

RXC*1.0 

Hl«(Z-0.)/Nl 

M«Nl/2 

SUM=RX0*RX1 

CCS5  J»1.M 

VVRR»0.0* (2*J-1 .  0)*Hl 

RX«OEXP(-(VVRR«*2)/2J 

SUH»SUM*4*RX 
85  CONTINUE 

IFIH-U  36  •36.115 
115  Mi»M-l 

CC65  J'l.Ml 

VVRR«0.0+2*J*H1 

RX*D£XPI-{ VVRR**2)/2) 

SUM»SUM*2*RX 
65  CONTINUE 
3t    TNsHl/3.*SUN/2.5066 

TNsO.S-TN 

GO  TO  ICO 
50  T.N-C5 

GO  TO  ICO 
IOC  WRITE(6.2)TN 

2  FORMAT!  ♦!'  .F10.51 
STOP 

ENC 
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APPENDIX  B 
Nonlinear  Models  Program  (NONLIN) 
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DEPARTMENT  CF  STATISTICS  AND  ThE  STATISTICAL  L.AECRATCRY  -  NCN-LINEAR  PROGRAM 

TITLE, CCST  FUNCTICN 
JC3  ..  XPRS0525 

DATA, 5. 

LABEL  (  11  "STRENGTH, STDS  ,R AG  lb S , A LPH AT  ,  CCST 

PARKLAE  1 1 )«Elf 82, E3» 84,35,86 , 8  7,  38 

METHOD, MARGUAROT 

ITERATIONS, 6C 

PARAMETERS, 8 

I N IT (1)«0. 7, 0.7, C.7,0. 7, 0.7, J. 7,0.7,0.1 

ObTPUT, STEPS 

MCDEL,5«Ffl,2,3,4] 

END  OF  PARAMETERS 

FCRMAT(F10.l,Fl0.2,FI0.1,F10.1,F10.l) 
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SU6RCUTINE  TRANS(X) 

ci.MEKSicr>   xii) 

RETURN 
END 


SLQRCITINE  FUNS(E,F.X) 

REAL«8  F,B(  1J 

REAL**  X( 100) 

COMMCN  NRAW  .NTRANSi  ICROPiNMCC 

F=(3(1)*X(1  )«*8  12)  )»(8I3)/(X(2)/X(1)  J**e(4J  )*(E(S)*M2]  **3(  6  J  )  *  (  B  { 
17)/X{4)**B (8) ) 

PLACE  FUNCTICNS  HERE 

RETURN 
ENO 


SU8RCITINE  LGRANG(8,GJ 

real*6  cm  ,ei i) 

.CCM-CN  NRAWfNTPANS,  ICPCP,  JvMCD 
COMMCN/CCNST/NITER,NPAR*rNCCNST,NVAR ,N,NCES 

PLACE  CCNSTPAINTS  HERE 


RETURN 
ENO 
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The  conventional    design  approach  based  on  a  "safety  factor"  gives 
little  indication  of  failure  probability  of  a  component.     In  reality,   the 
failure  probability  may  vary  from  a  low  to  an  intolerably  high  value  for 
the  same  safety  factor.     The  stress-strength  interference  theory  used  in 
this  thesis  is  a  better  method  to  calculate  the  reliability  of  a  com- 
ponent.    Time  dependent  stress-strength  interference  theory  is  also  dis- 
cussed. 

This   thesis   studies  an  optimization  problem  for  an  engineering  de- 
sign using  the  stress-strength  interference  theory.     The  problem  consists 
of  three  objectives,  which  are:     (1)  maximization  of  the  component  re- 
liability,  (2)  minimization  of  the  total   cost,  and  (3)  minimization  of 
the  weight.     Two  methods  of  multiple  objective  decision  making  are  used. 
They  are  nonlinear  goal   programming  and  the  Sequential   Multiobjective 
Problem  Solving  technique  (SEMOPS).     Goal    programming  requires  a  priori 
articulation  of  preference  information;  that  is,  the  decision-maker 
needs   to  give  goals  and  an  ordinal    ranking  of  objectives,  while  SEMOPS 
is  an  interactive  method. 

A  design  of  a  tensile  element  of  a  critical   component  of  a  system 
is   presented  as  an  illustrated  example.     The  coefficients  of  the  two 
objective  functions   (cost  and  weight)   in  this  example  are  estimated 
from  the  real   data  by  the  regression  technique.     The  optimization 
problem  is   formed  to  have   three  objectives  mentioned  before,  and  four 
decision  variables,  mean  tensile  strength,  standard  deviation  of  tensile 
strength,   radius  of  the  element,   and  manufacturing  tolerance  fraction. 
The  results  obtained  by  using  the  nonlinear  goal   programming  approach 


indicate  that  they  are  very  sensitive  to  the  ordinal   ranking  of  the  im- 
portance of  the  three  objectives  in  the  problem.     However,   the  results 
from  SEMOSP  are  different  from  nonlinear  goal   programming,  because  they 
indicate  that  they  come  from  the  trade-off  procedure  among  three  ob- 
jectives  in  an  interactive  manner. 


